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In the pyrochlore lattice Heisenberg antiferromagnet, for large spin length 5", the massive classical ground 
state degeneracy is partly lifted by the zero-point energy of quantum fluctuations at harmonic order in spin- 
waves. However, there remains an infinite manifold of degenerate collinear ground states, related by a gauge- 
like symmetry. We have extended the spin-wave calculation to quartic order, assuming a Gaussian variational 
wavefunction (equivalent to Hartree-Fock approximation). Quartic calculations do break the harmonic-order 
degeneracy of periodic ground states. The form of the effective Hamiltonian describing this splitting, which 
depends on loops, was fitted numerically and also rationalized analytically. We find a family of states that are 
still almost degenerate, being split by the term from loops of length 26. We also calculated the anharmonic 
terms for the checkerboard lattice, and discuss why it (as well as the kagome lattice) behave differently than the 
pyrochlore at anharmonic orders. 
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I. INTRODUCTION 

Highly frustrated magnetic systems are systems in which 
there is a zero temperature macroscopic classical ground state 

] degeneracyJi^ In experimental systems, this degeneracy is 
generically broken by secondary interactions, or by lattice 

] distortionsi^i^ However, even in toy models that include no 
such perturbations, one finds that the classical ground state 
degeneracy is broken by thermal fluctuations or quantum zero- 
point fluctuations. Such phenomena are collectively referred 
to as order by disorder^ 

Among three-dimensional geometrically frustrated sys- 
' tems, the most studied, by far, is the pyrochlore lattice, which 
is composed of the centers of the bonds of a diamond lattice, 
so the pyrochlore sites form corner sharing tetrahedra. Despite 
numerous studies designed to illuminate on the ground state 
properties of this model, in the large-S Iimi t^'^'^-^-'°-"''^''^ , a 
unique ground state has not been found for the pure, undis- 
torted, pyrochlore Heisenberg model. In this paper, we answer 
[ this question by finding the effective Hamiltonian that repre- 
sents the quantum zero-point energy to anharmonic order in 
spin waves. (A short report has appeared in Ref.fTil). 

We consider the nearest neighbor Heisenberg Hamiltonian 
on the pyrochlore lattice 

'■ 7^ = J^S, -Sj. (1.1) 

Here and below, (ij) denotes a sum over nearest neighbors. 
Classically, all states satisfying 

$^S, =0, (1.2) 

for all tetrahedra a are degenerate ground states, with energy 
— JNsS'^, where Ng is the number of spins (we reserve Greek 
indices for tetrahedra - diamond lattice sites- and roman in- 
dices for pyrochlore sites). 



A. Prior work 

In recent work,.2*l-L we have studied the quantum zero-point 
fluctuations of the large-S* limit of this model, and found that, 
to harmonic order in the 1/S expansion, there remains an infi- 
nite degeneracy of collinear spin states (although the entropy 
of this family is non-extensive). The degeneracy is associ- 
ated with an exact invariance of the harmonic order energy to 
a gaugelike transformation. Collinear configurations that are 
related by this symmetry have identical ^Mxei through all di- 
amond lattice loops, where the flux ipc through loop C with 
bond centers at (ii, Z2, . . . , i2n) is defined as 

The Ising variables rji = ±1 correspond to the classical spin 
direction along the collinearity axis. The harmonic ground 
states are all of the Ising configurations in one of these gauge 
families and we call them the -K-flux states, following Ref.Ua 
These consist of all collinear configurations whose fluxes 
through all hexagons (the shortest diamond-lattice loops) are 
negative: 

n (1-4) 

[The argument for ( |1.4t is given in Sec. |Vll] Some of these 
states are shown in Fig. 9 of Ref.fTl'. Furthermore, in Ref.fTll 
we constructed an effective Hamiltonian for the harmonic or- 
der zero-point energy, of the form 

Eltrn = NsS {Eo + X6$6 + A'8$8 + •••), d-S) 

where Eq, Kn are numerical coefficients that can be evaluated 
analyticallyii (Eq = -0.5640 , = 0.0136, Ks = -0.0033); 
here <I>2,i is the total flux (per lattice site) through all diamond- 
lattice loops of length 2n: 

^2n = ^ "^-c- (1-6) 

\C\=2u 
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In the interest of conciseness, throughout the rest of this 
paper we use the term state to mean "classical Ising configu- 
ration". In this paper, we go beyond the harmonic order in the 
expansion 1/S, to search for a unique semiclassical ground 
state, focusing in the asymptotic 5* ^ oo properties. We con- 
sider small quantum fluctuations about classical Ising config- 
urations such that the local collinear order is preserved. Our 
approach is aimed at deriving an effective Hamiltonian^ in 
terms of a much small number of degrees of freedom. 

Similar work has been previously done on the closely re- 
lated kagome lattice. This is a two-dimensional lattice, which 
is composed of corner sharing triangles. In the kagome 
Heisenberg antiferromagnet the zero-temperature classical 
ground states satisfy ( 11.21 ) for all triangles a, and harmonic or- 
der spin-wave fluctuations select all coplanar classical config- 
urations as degenerate ground states. A self-consistent anhar- 
monic theory breaks this degeneracy and selects one unique 
coplanar ground state -the so-called x state i'^''^-'^ 



B. Outline of the paper 

This paper is organized as follows: in Sec.|II]we derive the 
large-5 expansion by means of a Holstein-Primakoff trans- 
formation. We review some of the results of Ref. [HI on the 
harmonic theory. In Sec. |III1 we derive the mean-field Hamil- 
tonian for the anharmonic theory, and present a self-consistent 
variational scheme for solving it. 

Then, in Sec. |IV] we use a simple tractable example - the 
(tt, tt) state on the two-dimensional checkerboard lattice- in 
order to gain some analytic intuition on the behavior of the 
two-point correlation functions that governs the mean-field 
quartic energy, and the scaling laws involved. We find that 
these diverge as In S, resulting in anharmonic energy of or- 
der (InS*)^. In Sec. IIV Bl we argue that among all of the 
checkerboard lattice harmonic ground states, the quartic en- 
ergy is minimized in the (tt, tt) state, and show numerical re- 
sults to support this claim. We find that, due to the different 
symmetries of the checkerboard lattice and the Hamiltonian, 
the harmonic degeneracy in the checkerboard can be broken 
at the single-tetrahedron level, a result that cannot be carried 
over to the pyrochlore case. 

In Sec. |V] we present the main results of this paper - 
numerical calculations for the pyrochlore lattice. We find that, 
as in the checkerboard, the quartic energy scales as (In 5)^. 
We calculate the anharmonic energy for a large set of har- 
monic ground states and find that and that the anharmonic 
theory breaks the degeneracy between them. We derive ef- 
fective Hamiltonians for both the gauge-invariant and gauge- 
dependent terms in the quartic energy, and find a set of seem- 
ingly degenerate ground- states. 

Next, in Sec. |Vll we present a real-space loop expansion 
to explain the nature of the dominant term in the gauge- 
dependent effective Hamiltonian. We analytically derive an 
effective Hamiltanion, which is different from the one we con- 
jectured in the numerical fitting. Neverthless the leadng or- 
der terms of both effective Hamiltonians are minimized by 
the same set of states which, as far as we can tell, are all de- 



generate (both numerically and also to very high order in the 
effective Hamiltonian). 



II. SPIN- WAVE THEORY 

In this section, we expand the Hamiltonian jl.U in the 
semiclassical limit, in powers of 1/5. In Sees. IIIBjw e review 
some of the result in the harmonic theory of Ref. [Tlj, relevant 
to this paper. 



A. Large-5 expansion 

To study the quantum Heisenberg model, in the semiclas- 
sical limit of large S, we perform the Holstein Primakoff 
transformation. Since the harmonic ground states are all 
collinear,'' we shall in the following, limit ourselves to states 
in which each site is labeled by an Ising variable such that, 
without loss of generality, the classical spin is = i]iZ, and 
Siea ^ ^ ^'^^ ™^ tetrahedron a. Thus each tetrahedron in- 
cludes four satisfied - antiferromagnetic (AFM) - bonds and 
two unsatisfied - ferromagnetic (FM) - bonds. Notice that, 
whenever the spins satisfy the classical ground state condition 
( II. 2t . the sum of neighbor spins is (—2) times the spin on a 
site, i.e. 

= -2ry,. (2.1) 

jn.n. oH 

We first rotate the local coordinates to (?7ii, y, rjiz), and de- 
fine boson operators a^, a J such that 

St = VtiS-ala^), 

S+ = rnS'^ + iSy = ^j2S-ala, a,, 

S- = n,S^ -iSy = 4^j2S-ala,. (2.2) 

These operators satisfy the canonic bosonic commutation re- 
lations 

[ai,aj]=%, [ai,aj] = 0, [a|',a]]=0. (2.3) 

We now expand Eq. i2.2i in powers of 1/5, and express the 
Hamiltonian in terms of spin deviation operators 



< ^V^\|^ia^ + al), frf = -^^-(a,-a^), (2.4) 
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Wquart + 0(5"'), (2.5a) 
-JN,S'^ , (2.5b) 



-JSNs , 

J 
85^2 



(2.5c) 



where J = J(l + 1 /25). In the following, we shall set J —1. 
Somewhat redundantly, we also define Jij = J = 1 when 
are nearest neighbors, and zero otherwise (to simplify sums 
over just one index.) 

B. Harmonic Hamiltonian 

The use of the operators a^, allows us to represent the 
harmonic Hamiltonian (I2.5cb in block diagonal form 

Hha„„= (^^ ^) (^^^^-TrH, (2.6) 

where cr^ and cr*' are vector operators with respect to site in- 
dices, of length Ns, and the Ns x Ng matrix H has elements 

1 j 

Hij = \ \ I, j nearest neighbors . (2.7) 
otherwise 

The dependence on the particular classical ground state comes 
via the commutation relations 



(2.8) 



In Ref. luj we detailed the harmonic theory and the proper- 
ties of the eigenmodes. Here we briefly summarize the results 
relevant to this paper, for completeness. 



1. Diagonalization 

Define the Ng x Ng diagonal matrix n] by r]ij = rjiSij. 
Then spin-wave modes of any Hamiltonian of the form ( 12.61 ). 
with operator commutation relations ( 12. 8t are the eigenvectors 
{vm}, with eigenvalues {Am}, of the dynamical matrix rjfl: 



The eigenvectors satisfy a pseudo orthogonality constraint 

vjT]Vr,i(X 6lr,i. (2.10) 



The corresponding frequencies are tiuj„ 
the zero-point energy is 



= 5E(|Am|-l 



2S'|Am|, and thus 



(2.11) 



InRefs.OandU it was shown that the zero point energy is 
minimized for configurations that satisfy ( 11.41 ). A condensed 
version of this derivation shall be given later, in Sec. I VI Al 

For the Heisenberg Hamiltonian matrix ( 12.71 ) on the py- 
rochlore lattice, one finds that (for any Ising ground state) half 
the spin-wave modes have vanishing frequencies. These are 
the zero modes, which satisfy 



Vmii) = 0. 



(2.12) 



for all tetrahedra a. The two-point correlations (fluctuations) 
Gij of the spin deviation operators, it can be shown, are given 
by 



0. 



(2.13) 



It is clear from (12.131 ) that any mode v„i for which vJ„J7v,„ = 
0, exhibits divergent fluctuations. We call such a mode a di- 
vergent mode and it turns out that such a mode is necessarily 
a zero mode, i.e. A„i = 0. The converse is not true- most zero 
modes have nonsingular fluctuations. 



2. Ordinary modes 

The eigenmodes of Eq. i2.9i can be divided into two 
groups: half {Ns/2) of the modes have zero frequency. We 
call these generic zero modesr^ because the subspace that 
they span is identical for any collinear classical ground stated 
Since these modes have zero frequency, they do not contribute 
to the harmonic zero-point energy. 

The other half of the modes are called ordinary modesr^ 
and these modes can be naturally expressed in terms of 
diamond-lattice modes (recall that the diamond lattice has 
Ns/2 sites): an (un-normalized) ordinary mode v„i can be 
written down as 



(2.14) 



where the sum runs over the two tetrahedra to which site 
i belongs and is a vector of length Ns/2, living on the 
centers of tetrahedra (diamond lattice sites), and satisfying the 
spin-wave equation 

AmUm(a) = m(aP)Um{P) , (2.15) 

where the sum is over (diamond-lattice) nearest neighbors of 
a, and i{af3) is the pyrochlore site on the center of the bond 
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connecting a and /?. The diamond-lattice modes {u,„} are 
eigenmodes of an Hermitian matrix and therefore are orthog- 
onal to each other in the usual sense. We choose the normal- 
ization |u„i| = 1 without loss of generality. From (I2.14l i and 
( 12.15b one easily simplifies the pseudo-norm denominator in 
Eq. dTTSl l. 



(2.16) 



(valid only for ordinary modes). 

It is evident that the solutions of Eq. ( 12.151 1 are invariant 
under a gaugelike transformation of the state: if we trans- 
form ?7i —* TaTf^rji, where Ta,Tp = ±1, then the dispersion 
would not change, and the ordinary modes would transform 
'Um(a) TaU„i{a). Taking = —1 amount to flipping 
all of the spins in tetrahedron a. Such a transformation is not 
literally a gauge transformation since the flips must be corre- 
lated, so that the tetrahedron rule - J2i£a Vi — ^ from ( |1.2t - 
is preserved. Whenever two states are related by a gaugelike 
transformation, they have the same spin-wave eigenvalues 
and hence identical values of the total harmonic zero-point 
energy. 

Although most of the ordinary modes carry nonzero fre- 
quency, there is a subset of them that has A,„ = 0. It 
turns out that these are the divergent modes - modes that 
have vl^riVm = 0, and whose correlations are divergent [see 
Eqs. (Einil and ( I2l6] l]. 



3. Fourier transformed Hamiltonian 

In order to perform numerical calculations on large sys- 
tems, we must limit ourselves to periodic states. We shall 
assume a magnetic unit cell with A^j\/ sites arranged on a mag- 
netic lattice. In the simplest possible, Q = case, iVj\/ ~ 4 
and the magnetic lattice is the fee. Most of this work focuses 
on harmonic ground states, i.e., 7r-flux states. The smallest 
possible unit cell for that case has Nm = 16 sites. In practice, 
the calculation can often be simplified by utilizing the bond 
order, which may have a smaller unit cell3 

We Fourier transform the Hamiltonian i2.6i : 



H(q) 
H(q) 



TrH(q), 



(2.17) 



where cr^. 



cr^ are vectors of length Nm of the Fourier trans- 
formed X and y spin deviation operators. The wave vector q is 
in the Brillouin zone of the magnetic lattice. 




iq-(R+Ai) 



(2.18) 



where R is a magnetic lattice vector and I is a sublattice index, 
corresponding to a basis vector A;, i.e., for site i: = Ri + 



Upon diagonalization of the Hamiltonian [i.e., finding 
eigenmodes of r7H(q), where rj is now Nm x Nm], we ob- 
tain Nm bands within the Brillouin zone, half of which are 
of zero mode bands, and half are of ordinary modes. The di- 
vergent spin-wave modes can be shown to occur along lines 
in the Brillouin zone where an ordinary mode frequency goes 
to zero (we call these divergence lines) Each of these diver- 
gence lines is parallel to one of x, y, or z axes. 

The correlations of spin fluctuations can be expanded in 
terms of Fourier components, using Eq. ( |2.18t : 



N. 



M \ ^ 



N., 



q 



G;,;.(q)cosXi, • q, 



with 



(2.19) 



(2.20) 



where Zj and Ij are the sublattice indices of i and j, respec- 
tively, and Xjj =r,-rj - 



III. SELF-CONSISTENT ANHARMONIC THEORY 

This section develops our mean-field prescription to self- 
consistently calculate the anharmonic corrections to the en- 
ergy, for an arbitrary given state {rji}. First, (Sec. lIII A lb we 
decouple the quartic term Tiquart and write down a quadratic 
mean-field Hamiltonian. Next, we introduce a variational 
Hamiltonian as an approximation for mean-field problem 
(Sec, nil Bb . and in Sec. IIIIB~3] show that the variational form 
agrees with a general self-consistent approach in the large-5 
limit. In Sec. IIIICI we discuss how various fluctuations and 
energy scales depend on S. 



A. Decoupling scheme 

First let us work through the Hartree-Fock-like decoupling 
of the quartic term ( I2.5db of our spin-wave Hamiltonian^''. It 
turns out the decoupled coefficients depend on the (Ising) spin 
configuration in a simple fashion (Sec. lIII A 2b which allows 
us (in principle) to reduce the self-consistency conditions to a 
one-parameter equation. 



1. Energy expectation and decoupled Hamiltonian 

In a decoupling, one implicitly assumes a variational wave- 
function ^'mf, a priori unconstrained except for being Gaus- 
sian. Thus, it is specified by a harmonic effective Hamiltonian 
TCmf, defined so that 



{H 



harm 



n. 



quart , 



= {H 



mf; 



(3.1) 



where the expectations are taken with respect to Hmf itself. 

In light of Wick's theorem, we can immediately write the 
energy expectation by plugging into (|2.5cb and ( |2.5db the two- 
point correlations defined in (12.13b . but now using the Hmf 
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wavefunction: 

(Hharm) = 2 G,, + ^ - 5A^.) (3.2a) 

(3.2b) 



To make some expressions more compact, we define a bond 
variable. 



Tij = Gii — ijiTjjGi^ 



(3.3) 



Tij is, in general, «of symmetrioii and is defined only for 
(i, j) nearest neighbors (nonzero Jy ). 

Substituting ( I3.3l l into ( I3.2l i. and using ( 12.1b . we get 



(Hharm) = - ^1 [^''^J (^y + ^ 'ji) 



SNs 



Then 

-Emf = (TYmf) = ~ y~^(-^MF)ijGi 



(3.4a) 
(3.4b) 

(3.5) 



and [using ( 13.21 )1 we see indeed Hmf satisfies ( 13.11) . 

To write our decoupled Hamiltonian TYquait + ''^harm, we 
adopt a matrix form, in analogy with the harmonic Hamilto- 
nian ( 12.6b 



H 



MP 

H 





MF 



(3.6) 

defining the matrix elements in Eq. ( 13. 6b to depend on the 
correlations Gi,: 



{HMF)ij = 



-. Gii + Gjj — 2r]ir]jGi 



,(3.7a) 



(i/MF)„: = 1 + ^ X ^'^ ('mGj, - G^J ) . (3.7b) 



Recall from Sec. Ill Al that = 1 for nearest neighbors, oth- 
erwise zero. Thus, although G^ decays as a power law, Hmf 
has only on-site and nearest-neighbor terms. In terms of the 
Tij variables, eq. ( 13.71 ) reads 



{HMF)ij - 



(^^mf) 



(3.8a) 
(3.8b) 



All the machinery that was applied to H for the harmonic 
problem in Sec. IIIBI can now be applied to Hmf- In par- 
ticular, we can evaluate the correlations {Gy }, in terms of 
which the Hamiltonian matrix elements are written. Thus, 
by the self-consistent decoupling approximation we have re- 
placed the interacting spin-wave Hamiltonian by an effective 
non-interacting theory. 

Unfortunately, this does not yet give a solution, since the 
{Gij} are a priori unknown. We cannot just use the corre- 
lations obtained from the bare harmonic theory ( 12.61 ) for both 
practical reasons (G^ diverges in that case) and substantive 
ones: the theory would not be self-consistent - we would not 
recover the same correlations as those we put into it. A solu- 
tion may, in fact be obtained by successive iterations: assume 
a trial set of coefficients Hmf, compute the implied correla- 
tions, and define the next iteration of Hmf from ( I3.7ab . 



2. Simplified form of Tij and TimY 

In principle this iteration seems forbidding, but it is simpli- 
fied by an important fact, discovered numerically but verified 
analytically. For any Hmf approaching Tiharm, as should be 
the case for large S: 



AT 



I] ■ 



(3.9) 



Here F^^' and F^^' are diverging terms independent of i, j 
(and of the same order); whereas AF^ does depend on i and 
j. but is much smaller than F^^^ This was seen numerically 
in the outputs from a particular family of starting parame- 
ters, the family of variational wavefunctions ^'(e) specified 
by 'Hvai (e) [defined below in Sec. IIIIBll . More generally, an 
analytic explanation of the form ( 13.91 ), i.e. why F^ depends 
only on rji-qj at dominant order, is found in Appendix lAl [It 
follows from the gaugelike invariance, for the special case of 
Ising configurations that minimize the harmonic energy, the 
TT-flux states. One might crudely paraphrase that argument by 
saying the correlations that come out of the bare Hamiltonian 
have the form ([33 (albeit with divergent F^"', F^^)). 

Next, inserting the relation ( 13.91 ) into Eqs. ( 13.81 ), we can 
write the matrix elements of the mean-field Hamiltonian 



{HMF)ij 



52 



r(o) 



j2r('^?7«'?i3.10a) 



- ^(^r,, + AF,,) 



(i^MF). = (i-Jir(«)) + Ar(2) 



25 



(3.10b) 



To get the last line of Eq. ( I3.10bb . we used the z = 6 coordi- 
nation of the pyrochlore lattice, and the classical tetrahedron 
constraint X^iea Vi = ^ [from ( 11.21 )1. We now define 



J* 



52 



r(o) 



J* 



J Ji ' 



(3.11) 
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Note that I J* - 11 < 1. We obtain 



-^(Ar„- + Ar,o 



iHMF)ii ~ J* ( 1 + ^^out 



(3.12a) 



(3.12b) 



^2 



2S'2 



where 



^out - 



4r(2) 
52 J* 



(3.13) 



Thus, if we drop the much smaller terms in AF^ all the cor- 
rections are proportional to a single parameter F'^' times sim- 
ple functions of the spin configuration. 



B. Variational Hamiltonian 



The one-parameter dependence of Eq. ( |3.12t suggests we 
do not need to explore the full parameter space of trial Hamil- 
tonians to find the self-consistent mean-field Hamiltonian. 
Instead, we shall limit ourselves to a simplified variational 
Hamiltonian TYvar, which though it has just one variational 
parameter, appears to capture all the important properties of 
Hmf- (Specifically, Tivar approximates Hmf better and bet- 
ter in the limit S 00, as will be shown analytically below.) 

So, we wish to write a harmonic TYvar- as simple as pos- 
sible, to specify the Gaussian variational wavefunction 'I'var, 
its ground state (not necessarily equal to ^'mf)- Since Hmf 
- the solution to an unconstrained variational problem - has 
only nearest-neighbor terms, there is no loss of generality 
when we restrict our variational search to that form. [In 
contrasted, on the kagome lattice, the appropriate variational 
Hamiltonian had second- or third-nearest-neighbor (Heisen- 
berg) term a'^i'^-'^'^' , due to cubic terms in the spin-wave ex- 
pansion.] We thus adopt the simplest nontrivial form, the 
same as (13.6b . except with the diagonal block matrix Hmf 
replaced by 



Hvar = H + 6t]Ht] + el 



(3.14) 



where 6 and e are variational parameters. The 6 modifies the 
strength of AFM and FM bonds in opposite ways: namely, 
{H^ar)ij = (l+5)/2 for neighbors with 7]^ = and (i?var)ii = 
{l — 6)/2 for neighbors with r]i = This is the simplest 

possible form of a variational Hamiltonian that is consistent 
with the local spin symmetries. 

We do require invariance under global spin rotations, which 
means the Goldstone mode (associated with global rotation) 
must have zero energy. Its eigenvector vg has elements 



vaii) 



Vi. 



(3.15) 



Thus we require rjHvai vg = 0; inserting Eq. ( I3.7bl ) and writ- 
ing out each term, we first note Hvc = so our condition 
is 

= r]i Y,iH..rhvGU) - 4<5 + e (3.16) 



Thus (I3.14l i ends up having only one independent variational 
parameter e. It will become clear in the following, that the 
correct signs for the parameters are e > 0, 5 < 0. So, just 
writing out the components of Hvar as defined in ( I3.14l i. 
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(3.17a) 
(3.17b) 



Note that (ij) in ( |3.17a| i, and in similar equation pairs, applies 
only to nearest-neighbor sites. 

A more elaborate (multi-parameter) trial form of Hvar 
might improve the quality of the calculation, by exploring 
a larger set of variational wavefunctions; this is particularly 
important when the Ising configuration is not uniform from 
the gauge-invariant viewpoint (see Appendix IA 21 ). since ( 13. 9b 
breaks down in that case. Nevertheless, as we shall see nu- 
merically in Sec. [Vj the most important degeneracy-breaking 
effects are captured within this simple one-parameter theory. 



1. Self-consistent approach 



Revisiting eqs. ( 13.121 ), we see they reduce to Eqs. ( 13. 17b but 
with e Sout- Furthermore, as e ^ 0, it turns out £out(e) is 
increasing [indeed logarithmically divergent: see ( 14.171 ) and 
( |5.lb 1. So there is a unique self-consistent solution to 



£*sc = eout(e*sc) 



4F(2)(£*sc) 



(3.18) 



and at e 
get 



£*sc> [neglecting the AFy ) correction terms] we 



H 



MF 



J*Hv 



(3.19) 



Of course, the overall prefactor of J* has no effect on the 
spin correlations comprising Fy . Thus we have shown that, 
up to small corrections (of AF^ ), we in fact get out the same 
T^MF that we put in, so our theory is self-consistent. The only 
condition required for this to work was ( 13. 9b . 

Note F(") and F^^) are of order S In e, as will be explicitly 
verified analytically for the checkerboard lattice (Sec. IIV A 21 ) 
and the pyrochlore (Sec. IV Al l. The correction | AFy | in ( 13.91 ) 
is an order of magnitude smaller than F'^^' for all tractable 
values of e. 

If we had tried a different one-parameter form of variational 
Hamiltonian, where we add ±6 to the matrix elements H,, in 



a pattern other than the one in Eq. ( 13.141 ), the divergent Fy 
would indeed be regularized, but the dominant contribution 
would still be of the form ( 13.91 ). so self-consistency is lost: the 
output would not have the same as the input The only one- 
parameter nearest-neighbor variational Hamiltonian which is 
self-consistent is ( 13.14b . 
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2. Variational approach 

The above recipe is perfectly valid, but our actual calcu- 
lation was done somewhat differently. We diagonalized the 
Tivar to find a variational wavefunction '^^ar i^) and its corre- 
lations {Gij}, and computed an expectation Emf{£, S) [given 
by ( 13. 5H . We iteratively minimized Emf {^i S) with respect to 
e (for a given S), defining a unique optimal value e ~ e* (S). 

It will be shown below (in Sec. IIV A2I and IV Ab that 
e*{S) cx In 5/5. 



3. Equivalence of self-consistent and variational approaches 

It remains to be justified that e*sc^ defined self- 
consistently, should equal e*, defined by minimizing i?MF- 
This is expected, since the decoupling is variationally based: 
that is, a full variational optimization of Hmf with respect to 
all its parameters is equivalent to self-consistency with the de- 
coupling form, by construction. Thus, to the extent the full 
solution sticks within the subspace defined by Tivar (as we 
argued it did), the decoupling and variational minimization 
(both within that subspace) ought to agree with each other 

The test for whether our result really is self-consistent is 
that the diagonal elements ( |3.7b| i should be independent of i, 
and the off-diagonal elements (I3.7ab should depend solely on 
ijiTjj. Furthermore, we want {HMF)ij / {Hvai)ij to be equal for 
all i, j (for which Hij ^ 0). We indeed found (empirically) 
that this works when e = s*{S), i.e. [letting S*{e) be the 
inverse relation to £*(5)] 



variance^ 



(gMF(5*(£)), 
(iJvar(£))y 



(3.20) 



In Fig. [U we show an example of this for a particular state and 
a particular value of e. The crossing defines £*sc- in light of 
( 13.19b , but it is seen to happen exactly where £ = £*, thus 
empirically confirming the equivalence. 



C. Scaling 



Within the harmonic theory of Ch. IIIBI the fluctuations of 



the spin deviation operators scale as ((TjO-j 



0{S) 



omit the x and y component labels in these schematic expres- 
sions - and therefore we would naively expect, from the spin- 
wave expansion ( 12.51 ), that 



0{S) , (Hquart)r 



(3.21) 



However, Tiquart has an infinite expectation using the unmod- 
ified ground-state wavefunction of Tiharm, since the fluctua- 
tions diverge. Studies of the kagome lattic o'^''^''^ have taught 
us that, when anharmonic terms are treated self consistently, 
spin fluctuations of divergent modes are renormalized to finite 
values. In the kagome case {(Jidj ) = 0(5^/3) and the scaling 
relations are 




FIG. 1: (Color Online) Self-consistency of the matrix elements. We show 
the ratio of all nonzero elements of Hmf ™d Hvar (e) for the state shown in 
Fig. 9(d) of Ref.lTii Here e is set to 0.1. Each line represents a particular {ij) 
matrix element. Up to symmetiies of the configuration, there are 1 1 unique 
matrix elements for this state, some of which are virtually indistinguishable 
in the plot. All of the lines converge at 5*(e = 0.1) = 7.5 (up to a deviation 
which is much smaller than e). 



Note that the harmonic energy is not rescaled because the fre- 
quency of divergent zero modes is only 0(S^^^), which is 
negligible compared to non-zero modes' 0{S) frequency. 

One might expect the scaling ( 13. 22b to carry through to the 
pyrochlore lattice as welUi. However, the divergent modes 
of the kagome and the pyrochlore are rather different: in the 
kagome, due to the anisotropy between in-plane and out-of- 
plane spin fluctuations, all zero modes are divergent modes, 
so the kagome divergent modes span the entire Brillouin zone. 
In the pyrochlore, on the other hand, the divergent modes re- 
side only along lines in the Brillouin zone, hence the diver- 
gences (coming from these lines' vicinity) are weaker Be- 
low [see Eqs. ( 14.17b and (IS.lb l we shall find that this leads 
to logarithmic renormalization of the divergent fluctuations 



0{S In S), resulting in scaling 



{K 



quart/ 



MF 



(3.23) 



0(5) 



(Wquart)kag = 0(5^/3) . (3.22) be pOSitivC.] 



The singularity of the divergent modes' fluctuations, away 
from q = 0, is cut off by the variational parameter e. At 
q = 0, the divergence of (aidj) would be preserved, due to 
the physical Goldstone mode vq, but the Goldstone mode's 
contribution to vanishes such that the Goldstone mode 
does not contribute to the energy at any order in 1/5. 

Because it is technically difficult to deal with the diver- 
gence of Gij (q = 0) we shall, for now, retain both variational 
parameters. Thus we will have a handle on the fluctuations 
until we eventually take the limit 6 —e/A. [We find that 
0)^1/ V£ + 4(5, so that £ + 4(5 must be chosen to 
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IV. CHECKERBOARD LATTICE 

As a warm-up to the pyrochlore lattice problem, we 
first consider the same model on the closely related, two- 
dimensional checkerboard lattice. This case is more tractable, 
in that some expressions have a simple form which could not 
(or should not) be written out analytically in the pyrochlore 
case. 

The checkerboard lattice (see Fig. |2]) can be viewed as 
{001} projection of the pyrochlore lattice, and is often called 
the planar pyrochlore. The lattice structure is a square lat- 
tice with primitive vectors (1,1), (1,-1) and two sublattices 
corresponding to basis vectors (—1/2,0) and (1/2,0). We re- 
fer to the crossed squares as "tetrahedra" in analogy with the 
pyrochlore lattice, and we refer to any two sites within a tetra- 
hedron as "nearest neighbors" regardless of the actual bond 
length. 

Since the checkerboard lattice, as the pyrochlore, is com- 
posed of corner sharing tetrahedra, the derivation of Ch. 
remains valid. Note that we assume that all of the couplings 
within a tetrahedron are equal, even though in the checker- 
board lattice, the various bonds are not related by lattice sym- 
metries. Since the shortest loop in the checkerboard lattice is a 
square, the effective harmonic Hamiltonian for this lattice has 
the same form as the pyrochlore harmonic effective Hamilto- 
nian (11.5b . with the addition of a dominant term A'4$4, with 

^4 <0^ 

Thus, the harmonic ground states of the checkerboard lat- 
tice consist of all the zero-flux states, i.e., states with positive 
flux in all square plaquettes. Similar to the pyrochlore case, 
this is a family of states that are exactly degenerate to har- 
monic order, and in this case the residual entropy is 0{L), 
where L is the linear dimension of the systemjJi But since lat- 
tice does not respect the full symmetry of the tetrahedron, the 
selection effect of the anharmonic terms turns out quite dif- 
ferent (and essentially trivial) as compared to the pyrochlore 
case. 



A. The checkerboard (ttjTt) state 

One of the checkerboard harmonic ground states is sim- 
ple enough for the diagonalization of the variational Hamil- 
tonian ( 13.14b to be done analytically: the (tt, tt) state depicted 
in Fig.|2] In this state, the diagonal bonds in each tetrahedron 
are unsatisfied (FM), such that the symmetry of the lattice is 
conserved, and the magnetic unit cell has two sites. 



1. Harmonic Hamiltonian for checkerboard 

The Fourier transformed harmonic Hamiltonian for the 
(tt, tt) state is Eq. (12.17b . with 

/ 2cos2Q+ 2cosQ+cosQ_\ 

2cosg+cosg_ 2cos2g_ 



FIG. 2: (Color Online) The checkerboard lattice (tt, tt) state. The primitive 
vectors are the diagonal aiTows, and the piimitive unit cell is shown by the 
dashed square. The small arrows represents the two basis vectors. Here we 
show the (tt, tt) state: open (closed) circles denote up (down) spins. Dark 
(light) colored lines denote AFM (FM) bonds. 

where 

Q± = {q^±qy)/2. (4.2) 

The spin-wave modes can be found by diagonalizing the ma- 
trix J7H(q).— rj is a diagonal matrix with elements {r]i} along 
the diagonal (in our case 771 = 1, 772 = 2). Diagonalization of 
r)H.{c\) produces eigenmodes Vq and Uq for any wavevector 

q 

Vq = y^(cosQ+,-cosg_), Ay =/3q, 

= y5(cosg_,-cosg+), Ay -0, (4.3) 

satisfying the pseudo orthogonality condition V^rjUq = 0. 
The dispersions corresponding to Vq and Uq, respectively 
are 

Av, = /3q , Au, = , (4.4) 
Here we defined 

Qfq = 2(C0S^ g+ + COS^ g_) , 

/3q = 2(cos2g+ -cos2g_). (4.5) 

Thus, the ordinary spin-wave band has dispersion /kjq = 
2S'|/3q|, and the zero point energy can be easily calculated 

A^harm = ^ ^ /l^q - = N,S [ ^ - l] • (4.6) 

The fluctuations of the spin deviation operators (G'irn(q) = 
{a''^{l)a'-t^{m)) , where / and m are sublattice indices) can be 
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calculated from the spin-wave modes by Eq. ( 12.131 ) 



-7q 



-7q 



Here we used r||(q) [shown in Fig.[3ta)] for Fy (q), where 
both i and j are on the up-spin sublattice (and similarly for 
(A 7") ^ii- ^'^^ neighboring sites on different sublattices, we obtain 
[seeFig.Etb)] 



where 7q = 4cos(3_|_cos(3_, so that = y /3q + 7q- 
Eq. ( 14.7b shows that the fluctuations diverge wherever /3q van- 
ishes, i.e., along the lines in the Brillouin zone \Q+\ = \Q-\, 
which turn out to be g^; = or qy = 0. 



2. Anharmonic energy 

The variational Hamiltonian for the (tt, tt) checkerboard 
state is of the form ( 12.17b with the matrix ( 13.14b given by 

Hvar(q) = a^rnY^Ylr) + Sa^V^Vl + el , (4.8) 

Diagonalizing ?7Hvar(q), and keeping only the first order 
terms in 6, e results in Wq of order ^/e, along the diver- 
gence lines defined by /3q = 0, and a linear (in e,6) correction 
to ojq away from these lines. 

The fluctuations of the variational Hamiltonian are now: 



G(q) 



S 



aq(l + 5) + 2e -7q(l-(5) 

2D^{e,5)\ -7q(l-'5) aq(l+5) + 2£ 



Here we defined, for conciseness 



(4.9) 



Z?q(£, 5) = ^/32(1 - Sf + 4(aq + e)[a^5 + e) . (4.10) 

The fluctuations diverge (for nonzero e) only if /3q = and 
aq5+e = 0. If we take (5 ^ — e/4, to conserve the symmetries 
of the original Hamiltonian, we find one divergent mode: the 
q = Goldstone mode. 

In order to calculate the mean field energy ( 13. 5b . we are 
interested in combinations of the diagonal (on-site) and off- 
diagonal (nearest neighbor) fluctuations of the form F^ . We 
can write this as a sum over Fourier modes 



F, 



— E 

Nm ^ 



r.j(q) 



(4.11) 



[aq(l + 5) + 2e~ 7q(l - 5) cos q.^/y\ , 



(4.15) 

where F^^ (F|P is the bond variable for a bond oriented along 
the X (y) axis, connecting an up-spin and a down-spin. Note 
fliat Eqs. ( |4.13l ),( |4T4l ). and ( 14. 15b do not diverge at any value 
of q for £ + 4(5 = 0. Thus, we have regularized the fluctua- 
tions, and retained only one variational parameter. Since all 
sites are related by symmetry in this state, Fy = Tji. Fur- 
thermore F||(q) and F|j^(q) are related by a rotation of the 
Brillouin zone, and the real space correlations will be the same 
upon integration over the Brillouin zone. 

As we can see in Fig. [3] the divergent lines for F|| (q) and 
Fu (q) are both major axes, whereas F^^ (q) and F^j^ (q) only 
diverge along the y and x axes, respectively. Along the di- 
vergent lines, where (3^ = and ctq = |7q| = 4cos'^Q+, 
the values of the bond variables are, asymptotically Fy (q) = 
S\ sin 2Q-^\/2y/e. Away from the divergence fine. 



r,:,(q) 



51 sin2g^ 



2yFT4g2 



(4.16) 



where ^ 1 is transverse to the divergence line. Upon in- 
tegration of ( 14.131 ). ( 14.14b . and ( 14.15b over the Brillouin zone, 
the result is a logarithmic singularity in the fluctuations: 

45 

Ftt = = 2 ln£ + 0{e) = 2Ftx + 0{e) . (4.17) 

TT 

Observe that, in the notation of (iM . F^^ = Tn = F(°) + 
F(2) and F|| = F^o) - F^ ^), so the ratio 2 in Eq. ( |4l7] i is 
equivalent to the ratio 3 in ( IA12b . These fluctuations {Fy }, 
divergent as In e, enter quadratically into the anharmonic term 
of Eq. ( 13. 5b for the mean field energy Emf (The divergent part 
of the harmonic contribution, linear in {Fy }, cancels as was 
noted in Sec. 1111 A 21 ) 



with Tij (q) defined as 

Fy (q) EE G/.;, (q) - mVjGi.ij (q) cos ■ q . (4. 12) 

Here li, Ij are the sublattice indices of i and j, respectively, 
S^ij is the vector connecting the two sites. Nm is the number 
of points in the Brillouin zone, i.e. the number of sites in the 
magnetic lattice. 

In this case we obtain, for two neighboring sites on the same 
sublattice 



rTT(q) 

rii(q) 



5 



5 



[aq(l + (5) + 2£] sin2Q+,(4.13) 
[aq(l + <5) + 2£] sin2Q_.(4.14) 



Emf = £^harm + 5 X 0{s) - J2 (In £)' + 0(£ In £) 

= £;harm + 5x 0(£) + ^-^+0(£ln£). (4.18) 



Minimizing ( 14.181 ) with respect to e, for a given S^l (ig- 
noring the subdominant last term), we obtain e* (5) cx In 5/ 5 
and therefore the quartic energy Squart = Emf — -Bharm 
is quadratic in In 5. We remark that due to the logarithmic 
singularity, in a numerical calculation one would expect it 
to be hard to distinguish between terms of order C'((ln£)^), 
0(ln£), and 0(1) for numerically accessible values of e. 
Nevertheless, since we are doing a large-5 expansion, we are 
mostly interested in the asymptotic behavior 
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qx 71/ 2 



(a) 




qx 7T/2 



(b) 



FIG. 3: (Color Online) Bond variables in the Brillouin zone of the (tt, tt) 
checkerboard state, (a) Fij for two neighboring sites on the same sublattice. 
(b) Fij for two neighboring sites with riir]j = —1. In the case shown, the 
(ij) bond is along the x axis. The analytic forms of the functions are given 
in Eq. )4.13t and Eq. 14. 15t . respectively. 



B. Anharmonic ground state selection 



We start by explicitly finding the ordinary spin-wave modes 
of the harmonic Hamiltonian ( 12. 6t . Recall that the divergent 
modes are a subset (of measure zero) of the ordinary modes. 
Since we expect the divergent and nearly-divergent modes to 
dominate the fluctuations, we shall later limit ourselves to 
ordinary modes in the vicinity (in q-space) of the divergent 
modes. 

As we saw in Sec. Ill B 21 any ordinary mode can be 
written [Eq. (12. 14b 1 in terms of a vector u™, of length Ns/2, 
living on the centers of "tetrahedra". In the checkerboard case, 
these correspond to square lattice sites. {u,„} satisfy the spin- 
wave equation ( 12.15b . which can easily be solved by an ansatz 



(4.19) 



with j/q, = ±1 (to be determined). Plugging this into ( 12.151 ), 
we obtain, for any a 



^ ord 
Z ^ — ^ p n.n. of z 



iq-(r/3-ra) 



(4.20) 



As always, "ord" denotes a quantity limited to contribu- 
tions from ordinary modes. In order for the right-hand- side 
of (14.20b to be independent of a, we choose 



(4.21) 



It is easy to check that for (only) zero-flux states, the signs 
of \ya\ can be chosen consistently so that ( 14.21b is satisfied. 
(Note there is no need to assume the state is periodic.) Thus 
we obtain, from ( 14.20b . that for any checkerboard lattice zero- 
flux state, the dispersion is 



Aq = 2 COS qx COS qy . 



{A.22) 



Note that here and qy are shifted by (7r/2,7r/2) compared 
to Eq. ( 14.4b [for the (tt, tt) state]. This dispersion is shared by 
all of the harmonic ground states of the checkerboard. 

The (normalized) checkerboard-lattice ordinary spin-wave 
modes are, using ( 14. 19b in (12.14b . thus 



Now that we looked at the checkerboard (tt, tt) state, what 
can be said about the anharmonic selection in the checker- 
board lattice? The harmonic ground states in the checkerboard 
are the zero-flux state: all of the states that have a positive 
product over 77, around all square plaquettes. 

In this section, we shall first find the ordinary spin-wave 
modes (ignoring the generic zero modes, which are the same 
for all states), and then focus on the divergent modes to pre- 
dict which state is favored. Next, we show some numerical 
evidence to support are prediction. 



1. Spin-wave modes for a generic harmonic ground state 

In order to understand the leading order term in the anhar- 
monic energy, we restrict our discussion to the correlations 
due to divergent modes. We would like to derive an expres- 
sion for Tij, for any zero-flux state. 



Vae 



?q-rQ 



(4.23) 



The first term above vanishes upon summing over the lat- 
tice. 



2. Divergent correlations 



From (14. 23b . we can calculate the correlations due to ordi- 
nary modes, using (12.13b and ( 12.161 ) 

G°"^^J - n,iy ''o.vpacp (4.24) 



where 



_ S cosq • (r„ - r;3) 
q 



|Aq| 



(4.25) 
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is manifestly independent of which (zero-flux) state we have. 
Remember sum ( |4.24t has four terms; in the limit of a large 
system, the sum ( 14.251 ) converts to an integral in the stan- 
dard fashion. This is a special case of Appendix I A 1 1 ( 14.241 1 
corresponds to ( IA2b . and ( 14.25b corresponds to ( IA3b with 

gap = VaVpQaP- 

At this point it appears that we have a problem. The inte- 
grand in (14.25b diverges for any q along the divergence lines, 
and therefore we, of course, the correlations Ga, Gij diverge 
for the unperturbed harmonic theory. However, we have found 
that an adequate regularization scheme, such as the variational 
Hamiltonian (13.14b . cuts off the singularity and results in a 
logarithmic dependence. In particular, we have seen that, for 
the (tTjTt) state, X]q(VI'^ql) can be replaced by a constant 
C(e) which is logarithmic in e. Since the dispersion of Aq is 
the same for any zero-flux state, then C{e) can be assumed to 
be the same for all of the harmonic ground states. 

Without loss of generality, suppose site i is on the bond 
between diamond sites a and [3 and j is shared by a and /?'. 
Plugging this into Eq. ( 13.31 ) and using the relation ( 14.21b . we 
find the bond variables 

pord^^. ^ G°"d,,-,/,,;,G°^d^, (4.26) 

+7], cosq • {rp - Ta) - r/j cosq • (r^/ - r„)] . 

The last two terms in this expression are identically (since 
the sum is odd in q), and thus 

Assuming that the anharmonic selection is solely due to nearly 
divergent modes, we would like to focus on the vicinity of 
the divergence lines in the Brillouin zone: ~ ±7r/2 and 
Qy w ±7r/2. 

If the bond (ij) is diagonal, rp — rg/ = (±2, ±2), and the 
integral of the second term over any of the divergence lines 
is identically zeroi^ On the other hand, for a bond in the x 
(y) direction, the bond term in the bracket is +rjii]j for q = 
(±7r/2, Qy) (q — (q^;, ±7r/2)) and otherwise. 

Thus we find 

{SC{e) (ij) diagonal bond, 

5C(e)(l + ir;,77,) fe-)xorybond, (4.28) 
otherwise , 

Comparing to ( 13.9b . we see that r(°) = SC{e) while 
r^^) = ^r'"' on X or y bonds, but zero on diagonal bonds; 
the form is modified from ( 13.91 ) owing to the anisotropy of 
the "tetrahedron" in the checkerboard lattice (i.e., the inequiv- 
alence of the two kinds of bond.) 

Eq. ( 14.281 ) is by no means an exact result. We have made 
the following approximations in obtaining it: (i) Neglecting 
modes away from the divergence lines. This assumption is in- 
nocuous for large S, since the correlations are dominated by 
the vicinity of divergent modes. 




FIG. 4: (Color Online) The thi'ee possible polaiization axes for a single tetra- 
hedron. 



(ii) Neglecting all generic zero modes. In the checkerboard 
lattice, these modes, close to the divergence lines, can be 
shown to closely mimic the behavior of the ordinary modes, 
and will essentially increase C(e) by a factor of 2 (see Ap- 
pendix |A4l). 

(iii) Ignoring any additional effects due to the regularization 
scheme. Although this assumption is not a priori justified, 
we would like, as a first order approximation, to work with 
the bare harmonic Hamiltonian rather than the variational one, 
since it is easier to deal with analytically. We do not expect the 
regularization to qualitatively change the results we discussed 
in the following. 



3. Single tetrahedron 

To find the leading order quartic energy for a generic state, 
we consider the three possible bond configurations for a sin- 
gle tetrahedron, which can be viewed as three polarization 
axesi^isisl z (where all tetrahedra are oriented as in the (tt, tt) 
state), X and y (see Fig.|4]i;2^ 

Summing up the contributions, we obtain, for a single z 
polarized tetrahedron: 

<art = E ^^'^^r^^. « . (4.29) 

On the other hand, for x ox y polarization we find 

iSq^^^t « 2C{ef . (4.30) 

Note that in all cases ^ rjirjjV^^J^^ « to leading order, since 
the divergent modes do not contribute to the harmonic part of 
Emf in ( l33] l. 

Thus we found that the divergent contribution to the quartic 
energy is twice as large for x or y polarization as it is for z 
polarization. It follows that the effective Hamiltonian has the 
simplified form 

E^^,,, = Ns[A{S)-B(S)p.], (4.31) 

with B{S) w A{S)/2. Therefore the (7r,7r) state, in which 
all tetrahedra are z polarized, would be favored over all other 
zero-flux states, and thus is the unique ground state for the 
checkerboard lattice. 
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FIG. 5: (Color Online) A checkerboard lattice harmonic ground state. This 
state was constructed by flipping the bonds that cross each of the two dashed 
horizontal lines. 



2.5 




0.2 0.4 0.6 0.8 1 
Pz 



FIG. 6: (Color Onhne) Quartic energy for checkerboard lattice harmonic 
ground states. The energy -Equart is shown for e = 0.001, as a function 
of the fraction of z polarized tetrahedra, for various checkerboard lattice har- 
monic ground states. 



4. Numerics for full lattice 

To confirm Eq. ( 14.3 11 1 on the anharmonic selection among 
harmonic checkerboard ground states, we constructed various 
such states on a 8 x 8 cell (see Fig.|5]i in the following way: we 
started from the (tt, tt) state. There are 8 horizontal lines, that 
each go through the centers of 4 tetrahedra (dashed lines in 
Fig.|5]l. We choose any of the 2^ subsets of these 8 lines, and 
change the sign of rjirjj on every (vertical or diagonal) bond 
that crosses one of the chosen horizontal lines. It is easy to 
check, that each of these 2* transformations is a valid gauge- 
like transformation, since it does not violate the tetrahedron 
rule nor does it change the flux through any square plaque- 
tte. It turns out that of the 2^ that can be obtained, only 32 
are unique by lattice symmetry. Note that the construction of 
states, as well as our calculation, is based on bond-order^S 
and thus we need not worry about flipping an odd number of 
lines of this structurei^ See Ref. [Tl] for a detailed discussion 
of gaugelike transformations; for our purpose, it suffices to re- 
alize that each state that we generate is a valid classical ground 
state with zero flux through each plaquette. 

Whenever we flip a row of bonds, we change the polariza- 
tion of four tetrahedra from the z direction to the x direction. 
Based on the arguments of the previous section, we expect that 
the leading order term in the quartic energy would be propor- 
tional to the number of flipped rows. 

For each of these states, we calculate the quartic energy 
for a given value of e = 0.001, integrating over 41 x 41 
points in the Brillouin zone, equivalent to a system size of 
328 X 328, which is more than required to obtain good accu- 
racy (see Sec.[V]for more details about the numerical consid- 
erations). The results are presented in Fig.|6l as a function of 
the fraction of z-polarized tetrahedra pz- As expected we find: 
(i) the quartic energy is, for the most part, linear in pz- (ii) the 



energy span is of order 4(lne)^/7r''^. (iii) the ground state is 
the uniformly z polarized (tt, tt) state, (iv) the quartic energy 
of the (tt, tt) state is approximately half of the energy of the 
uniformly x polarized state. 

Given the clear differences in Squart {^iS) between the var- 
ious harmonic ground states, we expect that the same ordering 
would be conserved in the saddle point value i?quart {S) upon 
minimization with respect to e. Thus we can claim that the 
(tt, tt) state is the zero-temperature, large-S*, ground state of 
the checkerboard lattice model. This ground state is the same 
one found in large-A^ calculations for the large-S" limitJ2i2i 
The effective quartic Hamiltonian has the form ( 14.311 ) with 
the coefficients B{S) cx [In Sf and A{S) « 2B{S) to lead- 
ing order in S. We note that this effective Hamiltonian can be 
written in a more conventional form, in terms of Ising prod- 
ucts 

iS^Lt = N,A{S) - i?(^)^%.77, , (4.32) 

(y) 

where ^ ^ is a sum is over diagonal bonds only. 

The result is not very surprising: although we set the 
Heisenberg couplings to be the same for all bonds in the 
checkerboard lattice, there is no physical symmetry between 
the diagonal bonds and the non-diagonal bonds and therefore 
we should have expected to generate anharmonic terms con- 
sistent with the actual lattice symmetry. Thus, unfortunately, 
this does not provide a guide to lattices where all bonds in a 
tetrahedron are related by symmetry. 

V. EFFECTIVE HAMILTONIAN FOR THE PYROCHLORE 

We now turn our attention back to the pyrochlore lattice, 
where, due to the large sizes of the magnetic unit cells of 
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ground state candidates, it would be challenging, at the least, 
to do analytic calculations (as were done for the checkerboard 
in Sec. lIVI i. Since Sec. IIVI explicitly worked out the details, 
for that case, of implementing the self-consistent framework 
of Sec.|lII] we shall not belabor steps which are roughly paral- 
lel. However, the selection effects themselves — our ultimate 
motive — are quite different now, since the degeneracy is bro- 
ken by ordiinary modes in the checkerboard case. 

Our aim here is to calculate the quartic energy for a set of 
periodic states, and gather the energies we have calculated to 
construct an effective Hamiltonian. As seen in the harmonic 
theory of Ref. [Ill and in the large- theory of Ref. [l3, as 
well as the anisotropic perturbation theory of Refs. 125| and 
I26I it is natural that any non-trivial energy differences among 
states should be represented as a sum over loop operators. The 
effective Hamiltonian cannot take a local form: the local en- 
vironments that all spins see are the same. (Indeed, if we re- 
placed the diamond lattice by a (loop-free) 4-coordinated lat- 
tice Bethe, so that our spin sites formed a "Husimi cactus"—, 
then all Ising ground states would be equivalent by symme- 
tryii.) 

The numerical calculation is done as follows; for a given 
collinear classical ground state and a given value of e we di- 
agonalize the Fourier transform of the variational Hamilto- 
nian ( 13.141 ). keeping e + A5 infinitesimal. We find the bond 
variable Fy (q) for each wavevector on a grid of Brillouin 
zone points, and sum over these points to obtain F^ in real 
space. Once we have calculated Emf for many values of e 
(for a given collinear state), we can minimize it, for a given S, 
and find i?quart(S'). Our plan of action is to perform this nu- 
merical calculation of iJquart for a large database of collinear 
classical ground states and construct an effective Hamiltonian. 



A. Logarithmic divergences 

In performing the calculation, we find a distinct resem- 
blance to our findings on the checkerboard lattice: There are 
divergent modes along the x, y, and z axes in the Brillouin 
zone^ii and these modes dominate the mean field quartic en- 
ergy (and have no contribution to the harmonic order energy). 
The singularity of F;.;^ (q) is cut off, along the divergence 
lines, by a term of the order S/ %/e. The divergence peaks 
drop off to half of their maximum value at a (q) distance of 
order \/e, away from the divergence line. This means that the 
grid of wavevectors that we use must be denser in order to cap- 
ture the effect of the divergent modes, as e becomes smaller. 
Thus, we need to sum of the order of points, to obtain 
good accuracy. This limits the values of S that we can do the 
calculation for, and we have found no useful numerical tricks 
to get around it. Nevertheless, we can get results over about 
two orders of magnitude of S, which can be extrapolated to 
the 5 — > 00 limit. 

Upon numerical integration, we find, that as in the two- 
dimensional checkerboard lattice, the divergence of the fluc- 
tuations is logarithmic 

F„- oc |ln£|+0(e). (5.1) 



This numerical finding is somewhat surprising. We would 
naively expect that the bond variable Fy (q) would drop, away 
from the divergent lines, with a functional form ( |4.16t . as in 
the checkerboard. If so, as the transverse integration over q^ 
is now two-dimensional, the result would be a non-singular 
F 

It turns out that this expectation is incorrect because the dis- 
persion in the direction perpendicular to the divergence line is 
strongly anisotropic. For each value of q along the divergence 
line, there are two particular independent eigendirections of 
q^. For example, for a q = g-z divergence, the eigendirec- 
tions of qx are (1, 1, 0) and (1, —1, 0). If we call unit vec- 
tors along these eigendirections ei and §2, then we find that 
F„- cx 1/ ^£ + ((7_L • ei)2 + 1/ + {qi_ ■ G-^f . Integration 
over qx results in the logarithmic dependence on e of ( 15.11 ). 
as in the checkerboard case. In turn, as in Subsec. HV A 21 the 
logarithmic scaling of fluctuations in (15. Il l implies via (|3.18t 
that 

e*sc(5)(xi^^. (5.2) 

Finally, we know the decoupled quartic energy in Eq. ( 13.51 ) 
is a sum over products F^Fji, with Fy linear in InS*; since 
the divergent parts linear in F^ cancel out [as noted before 
dAl lb 1. the result is the anharmonic energy scales as (In 5)^, 
as announced in (13.231) . 



B. Gauge invariant terms 

For our database we calculated iJquai t on a sample of clas- 
sical ground states (not all of them 7r-flux states), that we 
constructed by hand, with unit cells ranging from 4 to 32 
sites. Two of these families consist of the zero-flux and 7r-flux 
states, which have uniform +1 and —1 products around all 
hexagons, respectively. In the other three gauge families, the 
hexagon fluxes are arranged in planes such that within each 
plane the flux is uniform; we call these the "OOOtt", "OttOtt", 
and "OOtttt" plane states, according to the stacking sequence. 

We minimize the Smf with respect to e at each value of 
S and obtain the energy shown in the inset of Fig. [7] We fo- 
cused on the five simplest gauge families. We minimize the 
i?MF with respect to e at each value of S and obtain the en- 
ergy shown in the inset of Fig. [T] We focused on the five 
simplest gauge families. We show the energies of all 16 dis- 
tinct Ising states belonging to the five gauge families. Due to 
the exact invariance of the (e = 0) harmonic energy under the 
gaugelike transformation, the total energies of states related 
by such transformations are, as expected, indistinguishable in 
the inset, since the harmonic term dominates. 

In the main part of Fig. [T] we show the anharmonic energy 
-Bquart for the Same states. As in the checkerboard lattice, 
the dominant part of the quartic energy is quadratic in In S, 
and of the order (InS")^. However, unlike the checkerboard 
lattice (compare to Fig.|6]l, we find that the energy differences 
between harmonically degenerate states are one to two orders 
of magnitude smaller than the dominant quartic energy. 
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We first consider the dominant gauge invariant contribution 
to the quartic energy. Since the invariants of the gaugelike 
transformation are products around loops, we search for an 
effective Hamiltonian in terms of the fluxes $2n, similar to 
the harmonic effective Hamiltonian (ll.Sl l. 



where we find, numerically 



10 



(5.3) 



Ao{S) « 0.300 + 0.0130(ln5)^ 

Ae{S) « -0.116 - 0.0030(ln5)^ 

^8(5) « -0.022 + 0.0055(ln 5)2, 

^10(5") w 0.008 - 0.0021(ln5)^ 



(5.4) 



Note that for large S, the signs of the coefficients Ae, As, and 
Aio are opposite to A'g, A'g, and A'lo in the harmonic Hamil- 
tonian. The differences in signs among the Ai{S) coefficients 
can explain why some of the lines in Fig. Qappear to be con- 
vex and other concave: each family of states is dominated by 
different flux loop lengths I. 

The gauge invariant terms can be heuristically explained in 
terms of the divergent modes: the quartic energy is large for 
states that have a large number of divergent modes. It turns 
ou t^'i^° that the number of divergent modes is linearly related 
to the flux terms $2?^: divergent modes proUferate to the ex- 
tent that the fluxes through loops of length 2n are (—1)". 

The above discussion of the gauge invariant quartic en- 
ergy ( I5.3l l is somewhat moot, inasmuch as it is negligible 
compared to the harmonic energy (|1.5l l, and it does not break 
the gaugelike symmetry. Nevertheless, one can clearly see in 
Fig.|2]that the anharmonic energy within each gauge family is 
not exactly the same, meaning that there is a gauge-dependent 
term in the variational anharmonic energy. 
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C. Gauge dependent terms and effective Hamiltonian 

Upon close inspection of Fig. [T] we see that some of the 
gauge families have a larger dispersion in their quartic ener- 
gies than others. But the quartic energy differences among 
the ground states of the harmonic Hamiltonian — the 7r-flux 
states — are much smaller than the gauge-invariant contribu- 
tion. We attribute this to the fact that, unlike the checker- 
board lattice harmonic ground states or even some pyrochlore 
gauge families, the 7r-flux states are completely uniform and 
isotropic (at the gauge-invariant level), and therefore there is 
no reason for the harmonic degeneracy to be broken at the 
single-tetrahedron level (see the discussion of Sec. IIV B 3l l. 
Indeed, in Appendix |A] we show that, the quartic energy due 
to ordinary modes of Tiharm - the dominant contribution - is 
gauge invariant among 7r-flux states. (This was not the case 
for the checkerboard case of Sec. IIV Bl ) We would expect any 
gauge-dependent terms in an effective Hamiltonian to not be 
as local as those in, say, Eq. (14.31b . 

In Fig. IS] we zoom in on the gauge dependent anharmonic 
energy, by showing the difference AE'quart = ^^quait - 
-E-quart, whcrc i?quart IS Calculated for 12 7r-flux states, and 



quart 



is the mean quartic energy of the states shown in the 



plot. 
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FIG. 8: Energy difference between i?quart of 12 harmonic ground states 
and the average of their energies Bquart- By taking differences between 
energies, we eliminate the (dominant) gauge-invariant term in the anharmonic 
energy. Each dashed line shows a fit in In 5^, for one of the states.. Note that 
there are several overlapping symbols along the bottom line, representing the 
degenerate states described later in the text (those with the maximum possible 
value of = Ns/3). 



FIG. 7: (Color Online) Q uai'tic energy i^quart for 16 classical coUineai' 
ground states. Equa,rt{S) was obtained in the variational calculation. The 
lines show a numerical quadratic fit in in S. Each gauge family (represented 
by 2-6 different states each) is denoted by a different symbol, of which tri- 
angles denote the harmonic ground states - the tt-Aux states. We show six 
vr-flux states, and their energies are virtually indistinguishable to the naked 
eye. The total energy i?MF is shown in the inset. 



In order to systematically search for a ground state con- 
figuration of the anharmonic effective Hamiltonian, we con- 
structed a large number of harmonic ground states using an al- 
gorithm for randomly generating gaugelike transformations jil 
Within unit cells that we used, of up to 192 sites, we believe 
that the algorithm performs an exhaustive search for harmonic 
ground states. About 350 states were found, inequivalent by 
lattice symmetries. (Notice that non-cubic cells were tried; 
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FIG. 9: (Color Online) The points with error bars are the numerical result 
Equa,Tt{S = 100) for 50 distinct tt-Hux states, which had been found us- 
ing our algorithm for generating gaugelike transformations. (Note that these 
energies are monotonic by construction, as the "state index" means simply 
the sequence when these energies are sorted. Every seventh energy is plot- 
ted. Shown for comparison are the energies predicted by the quartic effective 
Hamiltonian i5.5\ . using best-fit values for the three coefficients. The inset 
shows i?quart {S = 100) as a function of the effective Hamiltonian's leading 
term, Va- 



indeed, the optimal states described below require a cell di- 
mension that is a multiple of 3a/4 in the stacking direction.) 

The overall anharmonic energy (see Sec. Illl Cl l depends on 
S as (InS*)^, as does its gauge-invariant part [see Eq. ( 15.41 )1; 
is this also true for the gauge-dependent selection terms we 
seek? From what has been shown so far, that would be a plau- 
sible conjecture based on the scaling of the total energy, as 
well as the checkerboard case. Empirically, for each of our 
harmonic ground states, the S dependence of its energy (in- 
cluding the gauge dependent part) is well fitted by a linear 
or quadratic function In S (as seen in Fig. [8]). In fact, the 
checkerboard case is misleading: the anharmonic selection 
there (unlike the pyrochlore) depends on the ordinary spin- 
wave modes. The analytic derivation in Sec. |VT] shows the 
gauge-dependent term actually should scale as In 5*/ S; we do 
not understand the discrepancy between this and the numeri- 
cal results. 

In Fig. |9] we plot iJquait for the harmonic order ground 
states at = 100. There are two sources of error in this 
calculation; The first is the minimization error, represented 
by the error bars, which is due to the difference in energy be- 
tween consecutive value of e that we calculated, i.e. due to the 
"grid" in e-space. The second source of error is the grid used 
in integrating over the Brillouin zone, which is equivalent to a 
finite (albeit large) system size. This error becomes more sig- 
nificant for large values of S (i.e., smaller values of e), where 
the singularity of the divergence lines becomes narrower. The 
results shown are for 15'^ points in the Brillouin zone, for two 
different magnetic unit cells; a cubic 128 site unit cell, and a 
96 site tetragonal unit cell. 

As noted at the beginning of this section, we anticipate that 



an effective Hamiltonian should be represented by some sort 
of loop variables. We now consider an effective Hamiltonian 
of the form 

A£;^l^,t = Ce{S)Ve + CsiS)rs + Cio{S)Vio , (5.5) 

where Vi is equal to the number of loops of length I composed 
solely of satisfied AFM bonds. The form ( 15.5b was partly in- 
spired by the effective Hamiltonian from Ref. [l^ which is 
also a count of alternating loops (but with a broader definition 
of "loop" than here). Eq. ( 15.51 ) was guessed after fitting other 
forms with a variety of two- and four-spin terms involving the 
several closest neighbors. (Due to the ground-state constraint 
J2iea '^i ~ ^ ^'^'^ '■^^ TT-flux constraint ( 11.4b . there are numer- 
ous linear dependencies among such terms.) 

Also shown in Fig. |9] is a numerical fit to the effective 
Hamiltonian ( 15.51 ). For 5* = 100 we obtain 

Ce = -0.0621, 
Ca = -0.0223, 
Cio = -0.0046. 

(5.6) 

We ignore any constant terms here, as they belong in the 
gauge-invariant Hamiltonian ( 15.3b . 

While we cannot numerically repeat this calculation over 
a large range of values of S, in order to find the functional 
dependence C; (S) with good accuracy, we can obtain a rough 
fit by considering the small group of states depicted in Fig.[8] 
For these 12 states we obtain 

C'eiS) « -0.015 - 0.004(ln « 0.05 - 0.03 In S", 
C'siS) « 0.002 - 0.002(ln « 0.04 - 0.02 In S*, 
CioiS) « 0.0008 - 0.0005(lnS')2 sa 0.009-0.0041]^) 

Over our range of S" = 10 to 1000, either fit is plausible but 
In S" is a litttle better than (In S)'^. 

It must be noted that (at S ~ 100) the coefficients in ( 15.71 ) 
are bigger than ( 15.6b by nearly a factor of two; this is because 
the 12 states used were not sufficiently representative. Even 
though it is a rough fit, with significant error, it is clear (see 
the inset in Fig.|9]) that for a large number of states, the lead- 
ing order contribution to the energy is captured in Eq. ( 15.51 ). In 
particular, the numerical energy and the effective Hamiltonian 
agree as to which states have the minimum and maximum en- 
ergies. As it turns out, this can be predicted from the first term 
in ( 15.5b ; the highest energy states are those with the highest 
T'g value, namely Ns/&, which means 1/6 of all hexagons 
have alternating spin directions. It can be shown that, for tt- 
flux states, this is the smallest value that Vq can take3. The 
lowest energy states have Vq = Ns/i which is the highest 
possible value of Vq. 

D. Ground states 

Since the term is largest, and in view of the results just 
mentioned, it is a reasonable guess that the ground states are 
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(a) (b) 




FIG. 10: (Color Online) Projection of the slabs which form the near- 
degenerate anharmonic ground states of Ve, showing A layer in (a) and B 
layer in (b). The square shown is 2a X 2a. Open and filled circles represent 
spin up and down. Dashed lines are bonds outside the slab. In (b), one loop 
is outlined (numbered) from each of the two classes of hexagon mentioned in 
text and in Table H] there are also four classes for octagon placement, num- 
bered 3-6 here. 

a subset of the "hexagon-ground-states" that maximize just 
the Ve term. Since (see Appendix IC 1 b all hexagon-ground- 
states are degenerate at the octagon term too, only the much 
weaker 10-loop term might split these states, this assumption 
- confirmed numerically in the results of subsection IV CI - is 
very plausible. 

AU hexagon ground states found could be constructed by 
layering two-dimensional slabs (see Fig. [TOl i: they had unit 
cells of 48 spins (or larger). They were, within the numerical 
accuracy that we can obtain, degenerate for all values of S. In 
fact, we found these states share the same values of Vi for all 
loop lengths that we calculated (l < 16). Appendixlfijexplains 
these facts; indeed, it is shown that all loops are identical for 
/ < 26, and hence the stacked hexagon-ground-states must be 
exactly degenerate up to that order, at least for any effective 
Hamiltonian written in terms of loops [whether of the form 
i5.5i or the form to be derived in Sec. lVll . 

We conjecture that the stackings are, in fact, the only 
ground hexagon-ground-states, but this is unproven since we 
have not tried all possible unit cell shapes in the numerical 
enumeration. Appendix |C] explains how one could approach 
the ground state problem as a color-matching problem, but 
does not solve it. 

Although we shall find a different version of the effective 
Hamiltonian in Sec. IVH this section is valid for that too. All 
that matters is that the effective Hamiltonian depends on the 
Ising configurations of loops, and that the hexagon term dom- 
inates. 



VI. LOOP EXPANSION 

In Sec. |lll] we saw that in our self consistent theory, 
the mean-field Hamiltonian is proportional to the variational 
Hamiltonian 

HmF = J* Wvar (6.1) 

In fact, it turns out that the quartic selection effects of Hmf 
can be seen in the zero-point energy of Tivar, i.e. J* does 



not affect the selection. Therefore, we can try to understand 
the origin of the quartic effective Hamiltonian ( 15. St , by study- 
ing i^var, the zero point energy of the variational Hamilto- 
nian ( 13.141) . treated as a purely harmonic problem. 

In Refs. [Tl| and|^ we developed an effective Hamiltonian 
for the harmonic zero-point energy by a real-space loop ex- 
pansion. Below (Sec lVlST i. we shall use the same method as 
motivation for Eq. ( 15.5b . First, in Sec. ( IVI Ab . we shall give a 
quick summary of the results on Tiharm. Next, we represent 
the variational Hamiltonian in similar matrix notation, and re- 
peat the loop expansion (for the leading order in e), to derive 
an analytic effective Hamiltonian (Sec. IVICb . In Sec. IVI Dl 
we discuss the obtained effective Hamiltonian and compare it 
to the effective Hamiltonian we used in the numerical fit. 



A. Bare harmonic theory 

For this quick review of Ref. [U it will be convenient to 
rewrite some results of Sec. IIIBI using the matrix notation of 
( 12. 6b . as we note in each place. 

The spin-wave modes in the unperturbed harmonic theory 
are the eigenvectors of the equation [equivalent to (12.91 )1 

77Hv,„ = A,„v„ , (6.2) 

where H can be written as [equivalent to (I2.7b 1 

H=iwtW. (6.3) 

W is a Ns/2x Ns matrix whose (a, i) element is 1 if the 
pyrochlore site i is in tetrahedron a and zero otherwise. 

The spin-wave equation is transformed to the diamond lat- 
tice (which is easier to deal with, since it has fewer loops), by 
defining u„j = Wv,„. The diamond lattice modes satisfy the 
equation [equivalent to (12. 14b 1 

flUni = XmUm , (6.4) 

with the matrix //, = iWryW^. 

The elements of /j, only connect diamond-lattice nearest 
neighbors and are equal to the value of rj at the center of the 
bonds. 

_ / Viiap) a , /3 nearest neighbors , 
^"^5 ~ \ otherwise . 

As before, i{af3) is the pyrochlore site at the center of the 
diamond-bond (a/?). The zero point energy is 5 ^ |Am|, or 
in matrix notation 

Sharm=5Tr(i/x2)l/2_57V,. (6.6) 

For each a, the diagonal element (^/x^)q,q is equal to 1, and 
thus the square-root can formally be Taylor-expanded in pow- 
ers of fjL^ (or more exactly of /x^ — 41). 

£^harm/5' = 1 + ^ Q2„ Tvifl''^) - N, , (6.7) 
n=l 
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where the coefficients are 



the diagonal (with respect to the diamond-site index): 



Q2 



(-1)"- 



_i(2n-3)! 
8"n! 



(6.8) 



The details of the expansion were given in Ref. [Ill where 
the effective Hamiltonian ( 11.51 ). written in terms of {77;}, was 
derived from Eq. ( 16. 7t . However, the harmonic-order selec- 
tion can be explained with a "back-of-the-envelope" argu- 
ment, as in Ref.lgj: Tr /x^" is a sum of products of rjap over all 
closed paths in the diamond lattice. Since any path that goes 
back and forth is independent of {t?;}, the only paths that con- 
tribute non-trivial terms to the effective Hamiltonian are actual 
loops in the lattice. The first of these terms in for n = 3 (cor- 
responding to hexagons in the pyrochlore lattice). Thus, the 
first non-trivial term in the expansion favors states with neg- 



ative hexagon fluxes 
[Eq. (fL4l i]. 



the TT-flux states with rijeo = ^1 



B. Variational Hamiltonian 

The self-consistent theory (Sec. HiH i employs a variational 
Hamiltonian which has the same form as the harmonic Hamil- 
tonian but with H replaced by 



H^ar = H - ieryHjy + el 



(6.9) 



[to repeat ( |3.14| i and ( |3.l7| i1. Here e > is the (small) varia- 
tional parameter. The quartic energy is not equal to, but pro- 
portional to, the zero-point energy of the variational Hamil- 
tonian [with its parameter e* satisfying the self-consistency 
equation ( I3.l8b 1. Let us try to derive an expansion for this 
energy. 

The spin-wave modes are eigenvectors of the equation 



(6.10) 



Replacing H by ( 16.3b . we obtain 



A,„v,„ = Qr7WtW - ieWtWr? + £77 j v„ . (6.11) 

Clearly, the recipe for transposing this to the diamond lat- 
tice must be generalized to a more complex form than before 
(which must reduce to the old formulas in the case e = 0). 
Luckily, thanks to the simple form adopted for our variational 
Hamiltonian ( 13.141 ) it will suffice to expand the vector space 
of diamond modes from one to two components. Define the 
two vectors 



W?7v„ 



(6.12) 



For the case of e = 0, {u^} corresponds to ordinary modes 
and {u^„} to generic zero modes. 

It is convenient to introduce, analogous to /x, z/ = WW^; 
thus u is independent of {?/;} and has nonzero elements on 



4 a = /3 , 
v^p = i 1 a ,13 nearest neighbors , 
otherwise . 



(6.13) 



Still defining /x as in ( 16. 5b . we find [by multiplying 
Eq. ( 16.1 lb from the left by W and Wrj] the new equation 
of motion 



M 



with the 2Ns x 2Ns matrix M defined as 



M = 



(6.14) 



(6.15) 



The zero-point variational energy is 

E,,, = STv{hAy^^-SNs. (6.16) 

Note that now twice as many elements are summed in the trace 
as were in the bare harmonic version (16.61 ). One way to under- 
stand this is that the generic zero modes no longer have zero 
frequency and must explicitly appear in the zero-point sum 
5EIA™|. 



C. Expansion of variational energy 



The square root of ( 16.161 ) can be formally expanded in ex- 
actly the sum Eq. (16.71 ). but with the replacement /x^" 
M'^". In this trace expansion, each factor of /x or hops us to 
a neighboring site - with or without a factor of rjirjj, respec- 
tively - whereas a factor of 1 does nothing. We expect the 
lowest order non-trivial terms in the expansion to be of order 
6 in /X, v, since it takes (at least) that many hops to complete a 
hexagon, which is the smallest loop (in the pyrochlore lattice); 
these contributions come from the +Tr(M^) term 

Furthermore, since e is a small parameter, we shall expand 
the results in orders of e, keeping only the lowest order non- 
trivial term. Notice that for every 1 factor in ( 16.151 ). we pay 
the price of one power of e but do not gain a hop: hence, fac- 
tors of 11 cannot ever appear in a leading contribution. Such 
factors serve to "decorate" a basic loop, so that the same con- 
tribution reappears coming from higher powers of M and of 
higher order in e. They play a role similar to (and in addition 
to) the decorations by hops that retrace themselves, as found 
already in the bare harmonic theoryii. 

The upper-left block of M corresponds to {u^„} - the or- 
dinary modes, whereas the lower-right block corresponds to 
{u^} - generic zero modes (that acquire nonzero frequency 
in the variational Hamiltonian). Since the matrix elements of 
the sector always carry a factor e, the leading order terms 
in the small-e expansion will involve hops from the ordinary 
mode sector to the zero-mode sector and quickly return back. 
In this fashion, as conjectured in Appendix lAl we shall find 
explicitly that degeneracy breaking effects are due to the in- 
teraction between generic zero modes and ordinary modes. 
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All nonzero terms in a trace represent paths W of length 
21 < 2n on the diamond lattice that start and end on the same 
site (possibly retracing some bonds; also, 2n — 21 is the num- 
ber of factors Vaa which are diagonal with respect to sites. 
From here on we imagine having selected a particular path 
W, which can be expressed as a sequence of pyrochlore sites 
(diamond-lattice bonds) (ii, 12, . . . , i2;)' terms in the traces 
must be polynomials in the spins 77^^ , 77^3 , .... Then we consid- 
ering the terms due to Tr(M^") at each order in e. 

The leading order [0(1)] terms involve only the upper-left 
block (ordinary modes) of M. But it will be helpful to notice 
that Tr(/i^") = (2n)(/7w, where (pw = 11^=1 Vij > which gen- 
eralizes Eq. ( II.3I 1, to a general closed path. (Here the factor 2n 
accounts for different cyclic permutations inside the trace, i.e. 
different places the same loop could have been started. Note 
that any retraced portions in W have canceling contributions 
in the product </?w-) They are clearly gauge-invariant (See 
Appendix IA 2\ by the definition of the gauge-symmetry as de- 
scribed in Sec. II Al and are in fact exactly the same terms (/x^") 
that we had in the bare harmonic theory [Eq. ( 16.71 )1. Such 
terms in the effective Hamiltonian give the same value for 
all gauge-equivalent states, so they do not split the harmonic- 
order degeneracy and are not of interest here. 

In the next order, 0{e), we can have terms that take us out 
of the ordinary-mode sector in M and into the zero-mode sec- 
tor, but come immediately back. We obtain 

-^5ng2„eTr[/x2"-2(,y_8]l)«^] , (6.17) 

with the same 2n factor for cyclic permutations The trace in 
Eq. ( 16.171 ) contains two types of terms: Firstly, taking the site- 
diagonal {a = (3) element in each v, we obtain i(pw (where 
|>V| = 2n — 2.) As noted above, this is gauge-invariant hence 
not of interest. 

Secondly, taking the site-noti-diagonal elements of v, we 
obtain a products of all spins except two adjacent ones, i.e. 

fw'^VijVii+i , (6-18) 
3 

where we adopted the notation convention 111^^2^ = Vij ■ I" 
(only) the special case of a 7r-flux state, all products ipw along 
paths of the same topology are the same, and therefore a sum 
over all paths of length 2n amounts to a multiple of the clas- 
sical energy ViVj^ does not split any degeneracies. 
[More generally, within a family of non-7r-flux states, such 
terms do split the degeneracy and we must keep them. This 
is probably the reason that the dispersion of quartic energies 
among non-7r-flux states is notably larger than in the 7r-flux or 
0-flux states (see Fig.|7]i.] 

Moving on to the terms of order e^, we have contributions 
(i) from paths that hop once into the zero-mode sector (possi- 
bly staying there for at most one hop) (ii) paths that hop twice 
into the zero-mode sector (each time hopping back immedi- 



ately): 

^SnQ2ne'[TT [/x2"-3(zy -8l)fiu- 4.fi^''-^{i^ - 81)] 

2n-4 

+ - Tr - 81)«y/i2"-4-™(iy _ 8I)iv] | .(6.19) 

m=0 

The prefactor of 1 /2 in front of the second trace corrects the 
counting factor 2n, since each placement of the pair of vv> 
factors is counted twice in the sum. 

We now study Eq. (|6.19t . seeking to keep gauge dependent 
terms only. Start with the second term in the first trace, inside 
the curly brackets: — 4/i,^"^^(i/ — 81). In this term, only 
the site-diagonal elements in 1/ — 81 can contribute, since the 
path has to be of an even length. By the same arguments given 
above we just obtain (— 4)(4— 8)(y5vv which is gauge-invariant. 
Next, the first term in the first trace in Eq. ( 16.19b produces one 
gauge invariant term (for diagonal elements of v) plus one 
term that is gauge dependent: 

2n 

m Y VtkVik+2 = fwTw ■ (6.20) 
fc=i 

Every factor inside the trace involves a hop to a different site. 
Similarly, the sum over traces in (|6.19l l results, for a path W, 
in terms 

j=l k=j+2 

plus gauge-invariant terms that result from diagonal elements 
in 1/ — 41. This can be simplified into ^ipw{U^ — 2Ty^), 
where we define 

2n 

^^VikVik+i ■ (6-22) 

k=l 

Merging these two expressions together, we obtain, up to 
gauge invariant terms: 

-^SnQ2n£'^ ^ fwUw^ (6-23) 

|W|=2ri 

It is easy to see that only actual loops contribute interesting 
terms to Eq. ( 16.231 ) - all paths that go back and forth along the 
lattice add up to terms that are equal for all states that obey the 
"tetrahedron rule" ^^g^ = 0. Thus the anharmonic energy, 
to order e, can be expressed as a sum over lattice loops {£} 

i;var(gauge dep.) = "<52n Y "Pcpcl^ + 0{e^) . 

n=3 \C\=2n 

(6.24) 

Here, the coefficient Q2n is not quite the same as Q2n, since 
loop terms of length 2ji are renormalized by "decorated loops" 
of longer lengths. These are paths that go along the loop with 
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additional back-and-forth paths added to them. Such deco- 
rated loops have been discussed extensively, for related prob- 
lems, in Refs.flolfTll and can be summed up by use of simple 
combinatorics. 

Eq. (16.241 ) is the final result of this section and defines the 
quartic effective Hamiltonian E°^^j.^^. Assuming we chose 
e = s*{S), the self-consistent value, then each term in -Equart 
is oc Se"^, i.e. cx (In 5)^/5", in light of dS^ . We do not un- 
derstand the discrepancy (by a factor of 1/S) with with log- 
arithmic scaling of the fitted effective Hamiltonian in Fig. [8] 
and Eq. dSJb . 



to be the same for any of these states. Furthermore, at least 
for the stacked hexagon-ground-states found by the exhaustive 
search in Sec. [V] and described in Sec. erefsec:groundstates, 
many more terms are degenerate too. Each term appearing in 
Eq. ( 16.241 ) is the same in every state of this family, at least up 
to the terms for |£| = 16. Thus the degeneracy is broken only 
from a quite long loop that we anticipate to have a minuscule 
coefficient. 



VII. DISCUSSION 



D. Discussion of loop derivation 

With ( I6.24l i we can completely understand the essential 
features of the quartic effective Hamiltonian, and how the ana- 
lytic results of Sec.[yT]relate to the (prior) fit results of Sec.lVl 
Eqs. ( 16.241 ) and ( 15.51 ) are both sums over the same kinds of 
loops. The terms do not have the same analytic functional 
form, but are related, in being minimized by the same config- 
uration of alternating spins around that loop. Hence we under- 
stand how ( |6.24t and ( 15. 5t tend to be optimized by the same 
configurations, and hence why ( 15.51 ) was a good approxima- 
tion of the correct effective Hamiltonian. 

First, the leading order term in ( 16.241 ) is due to hexagons. 
Since the number of AFM bonds within a single hexagon (in 
a TT-flux state) can be 2, 4, or 6, and since \Uc\ = 2 is the 
same for both the case of 2 AFM bonds and the case of 4 
AFM bonds, then 

^^c\Uc\^ ^ -S2V6 + const . (6.25) 

O 

Thus, this term is in exact agreement the leading term in with 
Eq. i5.5i . It accounts for the largest contribution, sufficiently 
large that our ground state search can be limited to the subset 
optimizing the hexagon term minimizing | t/g P or equivalently 
maximizing Ve ■ 

The next to leading term is due to octagon loops. Already 
at this order, \Us\'^ is not independent of Vs- But, within tt- 
flux states, an octagon has (pc = +1, and since Qs = — l, then 
a large \Uc \ is favored. Clearly, a large X^Q l^-cl^ means a 
tendency to alternate and this correlates with large Vs, mean- 
ing that a large Vs is favored by Eq. ( 16.241 ). (In any case, 
among states optimizing ( 16.251 ), the octagon terms are always 
the same: see Appendix IC II ) 

As for loops of length 10 or longer, the situation is fur- 
ther complicated because the pyrochlore lattice has more than 
one kind (modulo symmetries) and Lpc may not be the same 
for different kinds of loop. Indeed, one kind of 10-loop has 
(pc = +1 while another kind has = — 1, in 7r-flux states. 
Therefore some of the 10-loops actually prefer to have a small 
\Uc\, and it is not certain a priori that Vw should be maxi- 
mized. 

But the role of larger loops simplifies in the special case 
of the hexagon-ground-states (the subset of 7r-flux states that 
optimizes Ve). The octagon terms (of either the fitted effec- 
tive Hamiltonian ( 15.5b or the analytic one ( 16.241 )) turn out 



We have calculated the anharmonic corrections to the spin- 
wave energy in the pyrochlore, and found that they break the 
degeneracy between the various harmonic ground states. We 
managed to numerically construct an effective Hamiltonian, 
and in Sec. I VII obtained an understanding of its terms. 

In retrospect, we should not have been surprised to find that 
the effective Hamiltonian is written in terms of loop variables. 
After all, in any collinear configuration, the local environment 
that each spin sees is the same for all sites. If the centers of the 
simplexes were put on a Bethe lattice rather than a diamond 
lattice, then all collinear configurations would be related by 
lattice symmetries and would therefore have the same energy 
(as was found explicitly in the harmonic theory of Re f. [Til 
and the large- theory of Ref.[l^ and in analogy to Ref.l27h. 
Thus any degeneracy-breaking terms must arise from lattice 
loops, so it is plausible that the effective Hamiltonian could 
be written explicitly in terms of loop configurations, but there 
are still multiple possibilities: the analytic derivation said the 
loop term is the square of the number of antiferromagnetic 
bonds along it [Eq. ( 16.241 )1 whereas a good numerical fit was 
obtained to a Hamiltonian that counts only the loops with all 
bonds antiferromagnetic [Eq. ( I5.5b 1. 

The anharmonic Hamiltonian is dominated by the smallest 
loops, the "hexagon" terms. The hexagon term's ground states 
are degenerate, having an 0{L) entropy; we conjectured that 
the stacked family in Sec. IV Di are all of its ground states, but 
we did not demonstrate it (see Appendix ICb. Within those 
states at least (and certainly to octagon order in any hexagon- 
ground-state), the count of many longer loops is constrained 
so that only a tiny term can break the degeneracy, which (for 
the stacked family at least) is only at the length 26 loops. To 
the accuracy layers of our numerics, all the stacked ground 
states are degenerate. 

What do our results say for realistic spins? First of all, 
the "small parameter" turned out to be l/ln5, which is not 
really small except at unphysical spin lengths [5* = 10-10"^ 
were used for numerical fits in Sec. IV CI . Still, our argument 
that only loop terms can break degeneracies still applies, so 
we expect the effective Hamiltonian takes similar functional 
forms for realistic S. It appears that only the first (hexagon) 
loop term will be important, since this will fix the values of 
the next few terms and only some very long loops will cause 
quite small splittings in these energies. So in practice this 
leaves a massive but non extensive degeneracy cxp(0(L)), as 
was already the case for the harmonic ground statsii (but with 
a smaller coefficient of L). 
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It is worth noting that the anharmonic selection effects in 
the pyrochlore turn out to be much weaker than in other 
closely related lattices: the two-dimensional checkerboard 
and kagome lattices. In the checkerboard lattice, which we 
discussed in Sec.|IV] many of the details are the same as in the 
pyrochlore: it is composed of corner sharing tetrahedra, the 
spin-wave Hamiltonian is the same, and the harmonic ground 
states are collinear states with uniform fluxes. Nevertheless, 
because of the anisotropy inherent to the two-dimensional 
checkerboard, the anharmonic energy breaks the harmonic de- 
generacy at the lowest order terms, of order (In S)'^. 

In the kagome lattice, the anharmonic selection is even 
stronger: first, there are cubic (in spin ct^/^) anharmonic 
spin-waves terms. In addition, because of the anisotropy be- 
tween in-plane and out-of -plane fluctuations about the copla- 
nar states, all harmonic zero modes possess divergent fluctu- 
ations and therefore the anharmonic energy scales as a power 
law in S-iSili^ 

Finally, we would like to mention that a similar calcula- 
tion can be carried out in the case of collinear states with 
nonzero magnetization, in the presence of a magnetic field. 
Such magnetization plateaus have been the subject of numer- 
ous recent studies . Our own harmonic work on the 
subject concluded that for a magnetic field that induces a 
collinear spin arrangement such that J^Vi = 2 in each tetra- 
hedron, the degenerate harmonic ground states are zero-flux 
states.— One could develop a self consistent variational treat- 
ment analogous to the one in this paper, to find that quartic 
ground state. Due to the asymmetry between f spins and j 
spins, there will be two independent variational parameters. 
In particular, the bond variables Fy are no longer expected to 
satisfy Eq. i3.9i . Rather, we expect the dominant terms in Fy 
to be F° + {rji - 77j )F(^) + rjir/jT^'^^ (see Appendix|A|. 



Acknowledgments 



1. Calculating correlations 

An ordinary mode Vm is a mode that can be expressed in 
terms of a diamond-lattice mode u„j by Eq. ( I2.14l i. The cor- 
relation function Gij was shown in Sec. Ill B Fl to be written as 
a sum over the spin-wave modes 

Yl 7:n :Vmii)vmij) ■ (ai) 

^ 2|vJ„?7v„| 

Restricting ourselves to the contribution of ordinary modes 
(denoted henceforth by superscript "ord"), and using 
Eqs. (IZlSl l and (|ZT6] |. 

= ViVj Sap- (A2) 

Q,/3:ea 

For (IA2b we defined, in analogy with (lAH 

^ — ^oid S 

9al3 = 2_^ 9U | "m(Q)Um(/^) ■ (A3) 

We need the bond variables (13. 3t , for a nearest-neighbor 
pair (ij), since that is how correlations enter our results [such 
as (I3.8l l1 . To express this for a particular pair, let a be the com- 
mon diamond site, and let (3 and /3' be the diamond sites at the 
far ends of the bonds on which sites i and j sit, respectively. 
Then 

5/3/3 + 9f}a - gp'a ^ 9pi3' ■ (A4) 

Note that the last line consists of one on-(diamond)-site corre- 
lation function, (the difference of) two nearest neighbor corre- 
lations, and one second-neighbor diamond mode correlation. 
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APPENDIX A: ORDINARY MODES 

To attempt to understand the results of the anharmonic cal- 
culation, the first thing we try is to calculate the contribution 
to the anharmonic energy due to ordinary modes, as we did, 
for the checkerboard lattice, in Sec. lIVBl The reason that we 
focus on ordinary modes is that, unlike generic zero-modes, 
we know how they transform under gaugelike transformation. 
In the checkerboard case, we saw (Sec. lIVBl ) that the anhar- 
monic selection can be explained in terms of the correlations 
due to ordinary modes in the harmonic Hamiltonian. As we 
shall see below, this is not true for the pyrochlore lattice, i.e. 
the ordinary modes produce a gauge-invariant quartic energy. 



2. Using the gaugelike symmetry 

Although we have been considering one particular classi- 
cal configuration, we can make use of the concept of gauge- 
like transformations (discussed in Sec. Ill B 2] ). The important 
points are the following: 

(i) Under a gaugelike transformation t (recall = ±1) 
the diamond-lattice spin-wave modes transform Um{ct) — > 

(ii) If two states have the same products of {rji} (flux) around 
each loop in the lattice, they are related by a gaugelike trans- 
formation. 

(iii) In particular, if the state has a uniform flux arrangement, 
(e.g. the TT-flux states), then any new configuration generated 
by a lattice-symmetry operation can alternatively be generated 
by a gaugelike transformation. 

The consequences of these points is that, for the 7r-flux states 
F'-^-' = is independent of a, (A5) 
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(since a gaugelike transformation would take a to (3 for any 
two diamond-sites a and /3). Similarly, it is easy to find that 
for nearest neighbor (diamond) sites a, (3 (sharing site i): 



r(i) 



Vi9at3, 



independent of i. 



(A6) 



and for next-nearest-neighbor (diamond) sites (3, (3' , con- 
nected by bond [ij): 



r(2) 



independent of (ij). 



(A7) 



In (IA7b . the sign was set so that F*^^' would be positive. Plug- 
ging these into (IA4b . we obtain 



(2) 



(A8) 



Since r°"^ must be invariant under a global spin-flip, we must 
have r'^) =0 and we obtain 



(A9) 



Eq. (|A9t is the key result of this appendix, the justification of 
Eq. It should be noted that T^"' andT^^' are both infinite 
in the bare harmonic theory, and are regularized by the varia- 
tional scheme. Here we assume that the regularization would 
not change the fact that F'"' and F'^^ are spatially invariant 
and gauge-independent. 

Furthermore, by the argument above, F'^''^ and F*^^^ are the 
same for any harmonic ground state (7r-flux state). Inserting 
Eq. iA9i into the mean-field energy ( 13. St , we quickly find that 
the ordinary modes' contribution to the anharmonic energy is 
gauge-invariant: 
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Note that the arguments above do not apply to the checker- 
board lattice, where all bonds are not equivalent by gauge- 
transformations - there is no transformation that can take a 
diagonal bond and turn it into a horizontal or vertical bond. 
Therefore, the correlations calculated from ordinary modes 
are sufficient to break the harmonic-order degeneracy in that 
case , as we find in Sec.HV] 



we could always do this well by using the wavefunction of the 
bare harmonic Tiharm- On the other hand, (13.41 ) [which is part 
of the expectation ( I3.5b 1 contains terms in {F^} which are 
divergent as e 0: these must cancel out, at the dominant 



order. In other words, F 



Tji, must cancel out. 



{H 



var /dominant 



^ rmi^^j + Tj.) « [r(") + f(2)] 

A^AFM[r(")-F(2)] (All) 



Since ( 13.91 ) says F^ (at dominant order) just depends on the 
sign of rjiTij , the sum groups into NpM terms for the FM bonds 
and A^AFM terms for the AFM bonds. But since A^afm = 
2Nym in any ground state. 



F(2)(£)/F(0)(£)^ 1 



(A12) 



valid for the limit e ^ 0. Numerically, F^2) appeared to be 
between F("V3 and F(o)/2, 



4. Role of generic zero modes 

Note that in the entire discussion, we have ignored the 
generic zero modes. Recall that divergent modes occur along 
lines in the Brillouin zone at q values for which the ordinary 
modes' frequency goes to zero. For q values close to these di- 
vergence lines, the zero-modes and small-frequency ordinary 
modes become close to each other (until they merge on the 
divergence lines; divergent modes are both ordinary and zero 
modes). The nearly divergent generic zero modes' contribu- 
tion to the correlations mirrors the contribution of the nearly 
divergent ordinary modes, and therefore F^ 
has the same functional form (IA9l l. 

In the self-consistent variational theory, the generic zero 
modes and the ordinary modes in the vicinity of the divergent 
lines interact strongly and, in fact, this interaction is responsi- 
ble for the degeneracy -breaking, as we observe in Sec.lVll 



2T°f and it 



APPENDIX B: STACKED GROUND STATES 

In this appendix, we analyze analytically the ground states 
of the effective Hamiltonians found in Sec. IV Cl and Sec. IVICI 
as summarized in Sec. IV Dl We assume a stacked spin config- 
uration (see Fig. [TOl i as this is what emerged from numerics; 
however, this is not yet proven. 



1. Layer stackings 



3. Relationof r'") toP'^' 

We take a moment to note that the parameters F^"-* and F'2) 
are not independent. We start from the variational Hamilto- 
nian (Sec. HITB] ). Notice that (Hvar) = -Eharm + 0{ e), O n 
the one hand, (TYvar) ~ ^-harm + 0(e), since [look at (13. 14b 1 



The pyrochlore sites can be broken into a stack of layers, 
each a/4 thick, where a is the lattice constant of the con- 
ventional cubic cell. The hexagon-ground-states are stack- 
ings of two kinds of slabs parallel to (say) the (001) plane: 
thin "A" slabs (thickness a/4) and thick "B" slabs (thickness 
a/2), which are stacked alternating A and B. A thin slab has 
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one level of chains along the [110] or ([10] direction, along 

which the spins repeat the pattern "H 1 — '. This pattern is 

reversed under a shift of [aOO] or [OaO], so the periodicity is 
\/2a X \/2a within a thin slab. 

A thick slab has two layers of spins, which form chains 
along the [110] and [110] directions, repeating the pattern 

"+ H ", such that the chain spins are parallel and the in- 

terlevel bonds are AFM in every tetrahedron spanning those 
two layers; within the thick slab, the spin pattern has a period 
2a X 2a. 

The inter-slab spin couplings cancel, so each slab has an in- 
dependent choice of two ways to align its spins. When there 
are m slabs of either kind, for a linear dimension in the stack- 
ing direction = m(3/4)a, the number of stacked spin 
states is thus 3 x 2^"' = 3 x 2^^^/^. This includes three pos- 
sible possible offsets (by multiples of a/4) in the z direction 
for the start of the stacking. [In a rectangular cell where L^. or 
Ly are also multiples of 3a (see below), we add similar terms 
counting possible spin stackings in the x oiy directions.] No- 
tice, apart that initial offset, the actual sites forming the layers 
are determined; only the spin directions are free. 

As a side remark, we can compare this to the family of har- 
monic ground states for the pyrochlore as described in Ref . § 
that was a stacking of only A slabs. The family of ground 
states of the effective Hamiltonian derived in the large-A^ the- 
ory for the pyrochlore'*^ is a stacking of alternating thin A' 
and B layers. The A' slab differs from the A layer shown in 
Fig. [Tot a) in that the spin patten is the same under a shift of 

Now we examine the slab stacking more carefully. The way 
a B layer adjoins A layers on opposite sides forces successive 
A layers to have opposite orientations: i.e., if one slab has 
chains along [110] the next one has them along [llG], etc. On 
the other hand, the way an A layer adjoins its neighboring B 
layers requires these B layers to have a relative shift in the 
xy plane of (a/4)[110] or (a/4)[liO] parallel to the A layer's 
chains. Hence, the xy offset of the B layer cycles through 
all four possible values in successive B slabs. The result is 
any periodic stack must have m even, e.g. m = 2 has a pe- 
riod [a/2, a/2, 3a/2] producing centered tetragonal cell. To 
directly repeat the same layer requires m to be a multiple of 
four, so the shortest cell (m = 4, = 3a) contains 12 layers 
of sites. 



2. Counting short loops 

Identifying ground states depends on counting the number 
of loops with various spin patterns, since this is what the ef- 
fective Hamiltonian depends on. We first do it for the shortest 
loops, starting with hexagons. A hexagon that satisfies the tt- 
flux constraint must have one of the four spin patterns shown 
in Table IJi; we label the types "i?2m" where 2n is the num- 
ber of AFM bonds in the loop. Also, independent of the spin 
pattern, the sites of a hexagon are placed in two possible ways 
within the layer stacking, which are the "classes" explained in 
the next paragraphs; the classes are also labeled in Fig.[TOl 

First, there are two classes of hexagon placement Class (1) 



TABLE I: Types of spin patterns in 7r-flux hexagon loops. Only 
hexagons with a loop product ipc = — 1 are included. Values are 
given for the two effective Hamiltonians, i5.5\ and ( 16.241 ) from the 
next section. 



Type class pattern 




H2 1 (+ + + ) 

m 1 (+ + + + +-) 
H4 2 (+ + - - +-) 
He 2 (+ - + - +-) 


4 
4 
4 
-1 36 



hexagons are centered on thin slabs. The two spins in the thin 
layer are opposite, and each pair within a thick layer is par- 
allel. Consequently, for each thin slab, the class 1 loops are 
half type H2 and half H2 (see Table Class (2) hexagons 
span one thick and one thin slab The part of the loop within 

the thick slab, always has H 1 — , so for each thick slab, the 

class 2 loops are half type H4 and half type Hq, of which the 
last is the type favored by the effective Hamiltonian. These are 
the four hexagon patterns satisfying the tt-Aux condition ( 11.4b ; 
that confirms that these slab stacked states are indeed har- 
monic ground states, a precondition for being hexagon ground 
states. Furthermore, since there are twice as many Class 2 
hexagons as Class 1, exactly 1/3 of all hexagons are type H2 
(the favored kind). Appendix F of ref. shows that a frac- 
tion 1/3 is the upper limit, so these are in fact hexagon ground 
states, too. A similar enumeration can be done of octagons. 
Again, for each particular type of spin pattern for an octagon, 
the number is the same for all our stacked hexagon ground 
states, therefore they are degenerate up to order 8. 



3. Long loops 

Symmetry can be used to show that much longer loops have 
the same count in all possible stackings. Say that a certain 
loop spans t slabs; the 2* possible spin states of those slabs 
are defined by (si, S2, st) where each Si = ±1 is a ref- 
erence spin in slab i. Now, a lattice symmetry operation g 
(which maps each layer to itself) has the action effect of flip- 
ping the spins in some slabs and not others: i.e. (si, S2, St) 
is multipled by some pattern of (71, 72, 7t) of ±1 factors, 
depending on g. Provided t is not too large, in fact every pos- 
sible pattern of 7^ is generated by some one of the lattice sym- 
metries: hence, all stacks of t slabs are related by symmetry 
and have the same counts of all possible loops. The smallest 
stack for which this no longer happens is when the first and 
last slab are stacked directly on top of each other, which (as 
worked out above) first happens for m = 4, meaning 1 2 layers 
or for i = 9 slabs (including the repeated one). The smallest 
loop which requires all of these slabs has length 2(12)-i-2 = 26. 

We conjecture that at order 26, the effective Hamiltonian 
does break the degeneracy. That will be a tiny energy: from 
( 15. 6t one could guess IC26 1 (for S = 100) is in the range 10"^ 
to 0{\Q~^) (depending whether one assumes an exponential 
decrease with 2n, or a power law). 
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TABLE 11: Supertetrahedra types: frequencies in hexagon-ground- 
states, and the counts of hexagons in each supertetrahedron (using the 
type labels of Table|I]) The types are given color names as explained 
in text. 



Type Name 


frequency 


Hexagons 






H2 , H2 ' 


H4 He 






orange 


white purple 


a purple 


1/3 





2 2 


b 








3 1 


c orange 


2/3 


2 


1 1 


d 





2 


2 



APPENDIX C: GROUND STATE PROBLEM AS COLORING 

Here we consider the ground states of the anharmonic ef- 
fective hexagon-order Hamiltonian, Ve- We review the ar- 
guments from Appendix F of Ref. The key idea is that, 
in a TT-flux state, there are constraints on spin arrangements 
due to the fact that different hexagons share edges. The 
level at which these contraints are first important is the super- 
tetrahedron, a cluster in the form of a truncated tetrahedron 
with four hexagonal faces. The centers of the super- tetrahedra 
form the complementary diamond lattice with the same lat- 
tice constant as the diamond lattice formed by centers of the 
original tetrahedron lattice. Each bond of the complementary 
diamond lattice (henceforth "superbonds") corresponds 1-to-l 
with a hexagon in the original pyrochlore lattice. 

We can classify supertetrahedra according to the types of 
hexagon loops appearing on their faces. Counting arguments 
there showed that there are four classes (Tablelllli and the total 
number of type 6 hexagons is maximized when only class (a) 
and (c) appear. 

1. Octagons in supertetrahedra 

First we can apply the supertetrahedron enumeration to 
show that all the hexagon-ground-states also are degenerate 
at the octagon term; we take advantage of the fact that ev- 
ery octagon is contained entirely within one supertetrahedron 
(three contained in each). 

We know that any hexagon ground state has fixed fractions 
of type (a) and type (b) supertetrahedra, as shown in Table HH 
But each of those supertetrahedra has a fixed pattern for its 



octagon loops: type (a) has one each of(H \ \ 1 1 — ), 

(+ + - H h-), and (+ H h - +-), while type (b) 

has one each of (+ H 1 1 — ), (+ + H h — ), and 

(+ + + + H 1 — ). Hence, any hexagon-ground-state has a 

fixed frequency of each octagon loop; from the list just given 
and the supertetrahedron frequencies in Table HH the octagon 
terms have the values Vs = 1/9, or mean [C/sP = 64/3. 
2. Node and superbond constraints as coloring rules 

A convenient necessary (though not sufficient) condition to 
be a hexagon ground state can be expressed as the following 
coloring problem on the complementary diamond lattice. For 
this purpose, the hexagon types (which are the superbonds on 
this lattice) are associated with colors, as are the supertetra- 
hedron types (nodes on the lattice). Then we have a com- 
plete covering by "purple trimers", consisting of two purple 
bonds (the middle node is purple and the other two nodes are 
orange. Simultaneously, we have a loop covering by orange 
loops (connecting orange nodes). Notice that, if we have such 
a coloring, we still must verify whether the can be filled in 
around each hexagon in a consistent fashion. 

In the stacking of Sec. IV Dl the supertetrahedra centered in 
B slabs are of type (a), and those centered between A and B 
slabs are of type (c). The purple trimer bonds are all oriented 
vertically (i.e. the three nodes are always at three different lev- 
els); this give 2^ degrees of freedom per B slab, accounting 
for all the spin entropy. The orange loops always run horizon- 
tally between the A and B slabs (perpendicular to the chains 
of that A slab). 

We conjectured, but did not prove, that the only hexagon 
ground states were the stackings of Sec. IV D| The special 
constraints of the stackings can be expressed, in the color lan- 
guage, as follows: 

(i) If a, (3, 7, 5 are four successive nodes connected by 
orange bonds, then the (a/?) and (7^) are oriented the 
same. 

(ii) If /3 is orange and 7 is a purple node, and {a(3) is 
the white bond into /3 while {j5) is the purple bond out 
of 7, then {af3) is never oriented the same as (7<5). 

We do not know if (i) and (ii) follow from the condition of 
having only type (a) and (c) super-tetrahedra, and so we do 
not know whether any hexagon ground state exists, besides 
the stacked family of Sec. IV Dl 
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niii] 



2|Aga;l sgn(cosgy) 



sm qy 

I cos I 



which, upon integration over qy vanishes for arbitrarily small 
^qx- 

3^ In Ref. [Ts] the polarization axis of checkerboard tetrahedra was 
denoted by a color Potts variable. 

3^ It should be noted, however, that when we impose a gauge- 
like transformation e.g. along a horizontal line, we are forced to 
change the vertical boundary condition from periodic to antiperi- 
odic or vice versa. 

3^ These states were illustrated in Fig. 3 of Ref. ]T(| with the stacking 
in the x direction, the B layers being those with all bonds AFM, 
or gray in that figure. 



